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Abstract 



Ansatzes for the Navier-Stokes field are described. These ansatzes reduce the Navier-Stokes 
equations to system of differential equations in three, two, and one independent variables. The 
large sets of exact solutions of the Navier-Stokes equations are constructed. 



AMS Mathematics Subject Classifications: 35Q30, 35A30, 35C05, 76D05 

O ■ 1 Introduction 

' The Navier-Stokes equations (NSEs) 

o ' 

which describe the motion of an incompressible viscous fluid are the basic equations of modern 
hydrodynamics. In ( |1 . 1| ) and below u = {u a (t, x)} denotes the velocity field of a fluid, p = p(t, x) 
denotes the pressure, x = {x a }, dt = d/dt, d a = d/dx a , V = {d a }, A = V • V is the Laplacian, 
the kinematic coefficient of viscosity and fluid density are set equal to unity. Repead indices denote 
summation whereby we consider the indices a, b to take on values in {1,2,3} and the indices i, j 
C~j ' to take on values in {1,2}. 



ut + (u- V)n - Au + Vp = 0, 

\ (1-1) 
divn = 



The problem of finding exact solutions of non- linear equations (LI) is an important but rather 
complicated one. There are some ways to solve it. Considerable progress in this field can be achieved 
by means of making use of a symmetry approach. Equations ( |1.1| ) have non-trivial symmetry 
properties. It was known long ago [|37| , |j that they are invariant under the eleven-parametric 
extended Galilei group. Let us denote it by Gi(l,3). This group includes the Galilei group and 
scale transformations. The Lie algebra ^4Gi(l,3) of Gi(l,3) is generated by the operators 

Po, Jabi D, P a , G a , 

where 

P = d t , D = 2td t + x a d a -u a d u *-2pd p , 
Jab = x a d b - x b d a + u a d u b - u b d u a , a ^ b, 
Ga = td a + d u a , P a = d a . 
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Relatively recently it was found by means of the Lie method || [§, |26|] that the maximal Lie 
invariance algebra (MIA) of the NSEs ([O]) is the infinite-dimensional algebra A(NS) with the 
basis elements 

d t , D, J ab , R(m), Z( X ), (1-2) 

where 

R(m) = R{m(t)) = m a (t)d a + m?(t)d u a - m? t (t)x a d p , (1.3) 

Z{x) = Z( X (t)) = X(t)d p , (1.4) 

m a = m a (t) and x = x{t) are arbitrary smooth functions of t (degree of their smoothness is 
discussed in Note |A.1|) . 

The algebra AG\(l,3) is a subalgebra of A(NS). Indeed, setting m a = S a b, where b is fixed, 
we obtain R(m) = db, and if m a = S a bt then R(fh) = G b . Here S a b is the Kronecker symbol 
( S a b = 1 if a = 6, Sab = if a^b). 

Operators ( |1.2| ) generate the following invariance transformations of system ( |1 . 1| ) : 

dt ■ u(t,x) = u[t + £,x), p(t, x) = p(t + e, x) 
(translations with respect tot), 

J a b : u(t,x) = Bu(t, B T x), p{t,x) = p(t, B T x) 
(space rotations), 

D : u(t, x) = e 6 u(e 2£ t, e 6 x), p(t, x) = e 2e p(e 2e t, e e x) 

(scale transformations), (1-5) 

R(m) : u(t,x) = u(t,x — rh(t)) + fht(t), 

p(t, x) = p(t, x — fn(t)) — rhu • x — \rh ■ rh t t 

(these transformations include the space translations 

and the Galilei transformations), 

Z(x): u(t,x) =u(t,x), p(t,x) =p(t,x) +x(t). 

Here e G E, B = {/3 a b} S 0(3), i.e. BB T = {S a b}, B T is the transposed matrix. 

Besides continuous transformations (|1.5|) the NSEs admit discrete transformations of the form 



t — t, X a — (2^6, Xb — Xbj 

h h ( 1-6 ) 

p = p, u a = u a , a ^ b, u = —u , 



where b is fixed. Invariance under transformations ( p..5| ) and ( |1.6[ ) means that (u,p) is a solution 
of ( |L~l| ) if (u, p) is a solution of (|T7T| ) . 

A complete review of exact solutions found for the NSEs before 1963 is contained in |l|]. We 
should like also to mark more modern reviews [|l(] : |?], |3(| despite their subjects slightly differ from 
subjects of our investigations. To find exact solutions of ( |1 . 1[ ) , symmetry approach in explicit form 
was used in j|, 31, 32,^, 20, 21,^, |I^,|l5, 12, 10, 11, ^] . This article is a continuation and a extention 
of our works [15, [l^, |l(], 11, 30]. In it we make symmetry reduction of the NSEs to systems of PDEs 
in three and two independent variables and to systems of ODEs, using subalgebraic structure of 
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A(NS). We investigate symmetry properties of the reduced systems of PDEs and construct exact 
solutions of the reduced systems of ODEs when it is possible. As a result, large classes of exact 
solutions of the NSEs are obtained. 

The reduction problem for the NSEs is to describe ansatzes of the form Q : 

u a = f ab (t, x)v b (u) + g a (t, x), p = f°(t, x)q(u) + g°(t, x) (1.7) 

that reduce system (1.1) in four independent variables to systems of differential equations in the 



functions v a and q depending on the variables uj = {w n } (n = 1,N), where N takes on a fixed 
value from the set {1,2,3}. In formulas ( |1.7[ ) f ab , g a , f°, g°, and oj n are smooth functions to 
be described. In such a general formulation the reduction problem is too complex to solve. But 
using Lie symmetry, some ansatzes (|1.7| ) reducing the NSEs can be obtained. According to the 
Lie method, first a complete set of A(A^5)-inequivalent subalgebras of dimension M = 4 — N is to 



be constructed. For N = 3, N = 2, and N = 1 such sets are given in Subsections A. 2, A.3, and 



A. 4, correspondingly. Knowing subalgebraic structure of A(NS), one can find explicit forms for 
the functions f ab , g a , f°, g°, and u n and obtain reduced systems in the functions v k and q. This 



is made in Sec. (N = 3), Sec. (N = 2) and Sec. H (N = 1). Moreover, in Subsec. 2.3 symmetry 



properties of all reduced systems of PDEs in three independent variables are investigated, and in 



Subsec. 4.3 exact solutions of the reduced systems of ODEs are constructed. Symmetry properties 
and exact solutions of some reduced systems of PDEs in two independent variables are discussed 
in Sections [5] and ^. In Sec. [?] we make symmetry reduction of a some reduced system of PDEs in 
three independent variables. 

In conclusion of the section, for convenience, we give some abbreviations, notations, and default 
rules used in this article. 

Abbreviations: 



the NSEs: the Navier-Stokes equations 

the MIA: the maximal Lie invariance algebra (of either a some equation or a some system of 
equations) 

a ODE: a ordinary differential equation 
a PDE: a partial differential equation 



Notations: 



C°°((to, ii), M): the set of infinite-differentiable functions from (to,t\) into K, where — oo < to < 
t\ < +oo 

C°°((tcb K. 3 )i the set of infinite-differentiable vector-functions from (to, t%) into M 3 , where — oo < 
to < h < +oo 

d t = d/d t , d a = d/d Xa , d u a = d/d u a, . . . 
Default rules: 

Repead indices denote summation whereby we consider the indices a, b to take on values in 
{1, 2, 3} and the indices i, j to take on values in {1,2}. 

All theorems on the MIAs of PDEs are proved by means of the standard Lie algorithm. 
Subscripts of functions denote differentiation. 
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2 Reduction of the Navier-Stokes equations to systems of PDEs 
in three independent variables 

2.1 Ansatzes of codimension one 

In this subsection we give ansatzes that reduce the NSEs to systems of PDEs in three independent 



variables. The ansatzes are constructed with the subalgebraic analysis of A(NS) ( see Subsec, A.2 
) by means of the method discribed in Sec.|B| . 

1. u 1 = |t| _1 / 2 (f 1 cost — v 2 sin r) + \x\t~ 1 — xx 2 t~ l , 
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\t\ 1 / 2 (v 1 sin r + v 2 cost) + \x 2 t 1 + xx\t 1 , 



u 3 = | t |-l/2 v 3 + ^ X3 t~ l , 

p = \t\~ x q + \>?t~ 2 r 2 + \r 2 x a x a , 



where 

yi = |t| _1 / 2 (xi cost + X2 sinr), y 2 = \t\~ l / 2 {—x\ sinr + x 2 cost), 

I/O 

2/3 = 1*1 x 3i K > 0, t = xln\t\. 
Here and below v a = v a (yi, y 2 , y 3 ), q = q(yuy^V3), r = {x\ + x 2 ,) 1 / 2 . 

2. u 1 = v 1 cos xt — v 2 sin xt — xx 2 , 
u 2 = v 1 sin xt + v 2 cos xt + xx\, 
u 3 = v 3 , 

p = q+l^r 2 , 

where 

y\ = x\ cos xt + x 2 sin xt, y 2 = —x\ sin xt + x 2 cos xt, 
y 3 =x 3 , x£{0;l}. 

3. u 1 = x\r~ 1 v 1 — x 2 r~ 1 v 2 + x\r~ 2 , 
u 2 = x 2 r~ 1 v 1 + xir~ 1 v 2 + x 2 r~ 2 , 

u 3 = v 3 + ri{t)r~ l v 2 + r)t(t) arctan x 2 /x\, 

p = q - 5r/ tt (t)(77(t)) _1 X3 - \r~ 2 + x(t) arctan x 2 /x±, 

where 

2/1 = *, m = r, y 3 = x 3 - rj(t) arctan x 2 /x\, 77,xeC°°((io,ii),M). 



(2.1) 



(2.2) 



(2.3) 



4 



Note 2.1 The expression for the pressure p from ansatz (2.S ) is indeterminate in the points t £ 



(to,ti) where n{t) = 0. // there are such points t, we will consider ansatz (2.S) on the intervals 
(tQ,ti) that are contained in the interval (to,t\) and that satisfy one of the conditions: 

a) n(t)^0 ViG(C*i); 

b) n(t)=0 Vt 

In the last case we consider iju/r] := 0. 

4. u = v t n l + (to • m)~ 1 v 3 m + (to • to) _1 (to • x)rht — yin\, 



q — |(m • m) {{fht ■ n l )yi) — (fh ■ rh) (fhu ■ x){m ■ x)+ 



(2.4) 



+ \{fhtt • r7i){fh ■ m) (in ■ x) , 



where 



yi = rt-x, y 3 = t, m,rf E C°°((t ,ti),R 3 ). 



n l ■ m = n 1 ■ n 2 = n\ ■ n 2 = 0, \n l 



1. 



(2.5) 



Note 2.2 There exist vector-functions n* which satisfy conditions ( \2.3[ ). They can be constructed 
in the following way: let us fix the vector- functions k 1 = k l (t) such that k % ■ fa = k 1 • k 2 = 0, 
\k l \ = 1, and set 



n 1 = k 1 cos ip(t) — k 2 s'mip(t) , 
n 2 = k 1 sin ip(t) + k 2 cos ip(t). 



(2.6) 



Then n] ■ n 2 = k\ ■ k 2 - ip t = if ip = J (kj ■ k 2 )dt. 



2.2 Reduced systems 



1-2. Substituting ansatzes (2^) and fl2,2| ) into the NSEs (|1 . 1| ) , we obtain reduced systems of PDEs 
with the same general form 



(2.7) 



Hereafter subscripts 1, 2, and 3 of functions denote differentiation with respect to yi, y2, and 2/3, 
accordingly. The constants ji take the values 



v a v\ 


- V L + H 


+ J1V 2 


= 


v a v 2 a 


~ V 2 aa + <?2 


- 7 1 t; 1 


= 


v a v\ 


- V L + 13 


= 0, 




< = 


72- 







1. 71 = -2x, 72 = -| if t> 0, 71 = 2k, 72 = § if t < 0. 

2. 71 = — 2x, 72 = 0. 
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(2. 



For ansatzes (2.3) and (|2.4j ) the reduced equations have the form 

3. v{ + v l v\ + v 3 vl - y 2 1 v 2 v 2 - {v\ 2 + (1 + n 2 y 2 2 )vl 3 )- 
-2ny 2 2 v 2 +q 2 = 0, 

v\ + v 1 vl + v^vl + y 2 l v l v 2 - [v\ 2 + (1 + v 2 V2 2 ) v h)+ 
+2r]y 2 2 vl + 2y 2 2 v 2 - rjy^q^ + %V% 1 = 0, 

v\ + v 1 v'l + v 3 v% - (uf 2 + (1 + v 2 y2 2 )4s) ~ 2 V 2 V2 3 4+ 
+2r ]1 y 2 1 v 2 + 2r ] y 2 1 (y 2 1 v 2 ) 2 + (1 + r] 2 y 2 2 )q 3 - 
-T} 11 r)- 1 y 3 - XW2 2 = °> 

V + wj + uf = 0. 

4. t>! + v j Vj - Vjj + qi + p l (y 3 )v 3 = 0, 

vl + t^'vj - v]j = 0, (2.9) 

vl + pHvs) = 0, 

where 

/°* = P i (V3) = 2(rh ■ m)~ 1 (rht ■ ft 1 ), 
p3 = p 3 (y 3 ) = (m • m)~ 1 (rht • m). 



(2.10) 



2.3 Symmetry of reduced systems 

Let us study symmetry properties of systems (p77|) , (|2.8[) , and (^fi|) . All results of this subsection 



are obtained by means of the standard Lie algorithm |28|, 27]. First, let us consider system (2/7). 

Theorem 2.1 The MIA of system (\2. % ) is the algebra 

a) < d a ,d q ,j{ 2 > if 71 / 0; 

b) <d a ,d q , Jl b > if 7l =0, 72 /0; 

C ) <d a ,d q ,4 b ,Dl > if 7l = 72 = o. 

Here J\ b = y a d h - y b d a + v a d v b - v b d v a , 
D\ = y a d a - v a d v a - 2qd q . 



Note 2.3 All Lie symmetry operators of ( \2. T ) are induced by operators from A(NS): The operators 
J^b and D\ are induced by J a b and D. The operators c a d a (c a = const) and d q are induced by 
either 

R{\t\ l l 2 (c\ cos t — C2 sinr, c\ sinr + C2 cos r, C3)), 



where r = xln|t|, for ansats (2.1) or 

R{c\ cos xt — C2 sin xt, c\ sinxtf + c 2 cos xt, c 3 ), Z(l) 

for ansatz ( \2.Q ), respectively. Therefore, Lie reductions of system \ ) give only solutions that can 
be obtained by reducing the NSEs with two- and three-dimensional subalgebras of A(NS). 
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Let us continue to system (|2.8| ). We denote A max as the MIA of (2.8). Studying symmetry 
properties of ( |2.8[ ), one has to consider the following cases: 
A- r],x = 0. Then 

A^= <d \ D lR 1 ^y 1 )lz 1 (X(y 1 )) >, 

where D\ = 2yid\ + y 2 d 2 + y^dj, - v a d v a - 2qd q , 

#1(^(2/1)) = ^3 + tpid v s - Tpimdq, Z 1 (X(y 1 )) = X(yi)d q . 

Here and below if) = ip(yi) and A = X(yi) are arbitrary smooth functions of y± = t. 

B. 77 = 0, x 0- I n this case an extension of A max exists for \ = (Ciyi + C 2 )~ l , where 
C\,C 2 = const. Let C\ 7^ 0. We can make C 2 vanish by means of equivalence transformation 
( A. 6 ), i.e., x = Cyf 1 ' where C = const. Then 

A™*=< D l 2 ,R 1 (^y 1 )),Z\X(y 1 ))>. 
If C\ = 0, x = C = const and 

A max =< d l ,R l ^(y l )),Z\X{y l ))> . 
For other values of x, i.e., when X11X 7^ XlXi, 

A ma X=< fl 1 (^(y 1 )),Z 1 (A(yi))>. 



C. 7/ 7^ 0. By means of equivalence transformation ( |A.6|) we make x = 0. In this case an 
extension of A max exists for rj = ±|Ciyi + C2I 1 / 2 , where C\,C 2 = const. Let C\ 7^ 0. We can make 
C 2 vanish by means of equivalence transformation ( |A.6| ), i.e., 77 = C|yi| 1//2 , where C = const. Then 

A ma X=< D l } R2 Q yi \l/2 );R2 Q yi \l/2 l n \ yi \) , Z 1 (X( Vl )) > , 

where R 2 (ip(yi)) = tpd3 + tpid v 3. If C\ = 0, i.e., r\ = C = const, 

A max =< d 1 ,d 3 ,y 1 d 3 + d v3 Z\X(y l )) > . 
For other values of rj, i.e., when (ry 2 )n 7^ 0, 

A ma X =< i?2 ( 7? ( yi )) ) i ?2 ( r? (y 1 ) f( v (y 1 ))- 2 dy 1 ),Z 1 (X(y 1 )) > . 



Note 2.4 In all cases considered above the Lie symmetry operators of §2Jy are induced by operators 
from A(NS): The operators d\, D\, and Z 1 (X(y\)) are induced by dt, D, and Z(X(t)), respectively. 
The operator R(0, 0, ip(t)) induces the operator i?i (-0(2/1 )) for r] = and the operator R 2 (ip(yi)) (if 



tpurj — V"7n = 0) for rj 7^ 0. Therefore, the Lie reduction of system j\2.^) gives only solutions that 



can be obtained by reducing the NSEs with two- and three- dimentional subalgebras of A(NS). 



When 77 = x = 0, system ( |2.8j ) describes axially symmetric motion of a fluid and can be 
transformed into a system of two equations for a stream function vl/ 1 and a function vf 2 that are 
determined by 

*3 = V2v\ *2 = -V2V 3 , = y 2 v 2 . 



The transformed system was studied by L.V. Kapitanskiy [20, ^jf| . 
Consider system (|2.9[) . Let us introduce the notations 



* = 2/3, p = p(t) = f p 3 (t)dt, 
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z\\(t)) = \(t)d q , s = d v3 - P l {t) yi d q , 

E (x(t)) = 2xd t + XtVidy, + (xttVi ~ Xtv l )d v i - (2x t q + \xtttVjyj)d q , 
J\2 = V\di ~ Vidi + v l d v 2 - v 2 d v i. 



Theorem 2.2 The MIA offing) is the algebra 

1) < R 3 (^(t), *P 2 (t)), Z l (X(t)), S, E(xHt)), E( X 2 (t)), v 3 d v3 , J\ 2 >, 
where x 1 = e~ p(t) J e p ^dt and X 2 = e~ p{t \ if p i = 0; 

2) < Rsi^it), V 2 (0), Z 1 (X(t)), S, E( X (t)) + 2 ai v 3 d v3 + 2a 2 j\ 2 >, where 
a\, a2, and 03 are fixed constants, x = e~ p ^ f e p ^dt + 03^ , if 

1 - e %Pp-\- a ^ (Ci cos(a 2 In p) - C 2 sin(a 2 In p)) , 



P 



p 2 = e 2 p p 2 ai ( Ci sin(a 2 In /5) + C 2 cos(a 2 hip) 



with p = p(t) = \Je p ^dt + a 3 \, Ci,C 2 = const, (Ci, C 2 ) / (0, 0); 

3) < R^it), tp 2 (t)), Z 1 (X(t)), S, E( X (t)) + 2a 1 v 3 d v a + 2a 2 j\ 2 >, where 
at and a 2 are fixed constants, x = e~ p ^ , if 

p i = e lp-*ip( Cl cos(a 2 p) - C 2 sin(a 2 p)) , 
p 2 = el p - aip (d sin(a 2 p) + C 2 cos(a 2/ 5)) 

with p = p(t) = J e p ^dt, Ci,C 2 = const, (d, C 2 ) / (0, 0); 

4) < R^it), i) 2 {t)), Z 1 (X(t)), S > in all other cases. 

Here ip % = tp % (t), A = X(t) are arbitrary smooth function of t = y%. 



Note 2.5 If functions p b are determined by t \2.1L ), then e p ^ = C\m(t)\, where C = const, and 



the condition p % = implies that m = \m(t)\e, where e = const and |e| = 1. 



Note 2.6 The vector-functions n l from Note 2A are determined up to the transformation 

n 1 = it 1 cos 5 — n 2 sin 5, n 2 = n 1 sin 5 + n 2 cos <5, 
where 5 = const. Therefore, 5 can be chosen such that C 2 = (then C\ 7^ 0). 
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Note 2.7 The operators Rz(ip ,ip 2 ) + aS and Z 1 (X) are induced by R(l) + Z(x) and Z(X), respec- 
tively. Here I = Tp l n l + ip 3 rh, ip 3 (fh ■ rh) + 2ip % (n^. • rh) = a, 

X — i(m • fh)~ l {(rht ■ h 1 )^ 1 ) 2 — \{m t t • h t )ip 3 ip l + \{ht • h % )^ % = 0. 

If rh = \rh\e, where e = const and |e| = 1, the operator j\ 2 ^ s induced by e 1 J23 + e 2 J31 + e 3 Ji 2 . 
For 

rh = f3 3 e at (/3 2 cos r, f3 2 sin t,/3i) t 



with t = xt + 5 and (3 a = const, where ft 2 + fi\ = 1j ^ e operator dt + xJi 2 induces the operator 
y 3 - PlxJl2 



d V3 — 0iftJi2 + crv 3 d v i if the following vector- functions n % are chosen: 



h 1 = k 1 cos f3\T + k smf3\T, h 2 = — k sin/?ir + k cos/?it, (2-11) 

where k 1 = (— sinr, cos r, 0) T and k 2 = cost, /?i sin r, — /3 2 ) T . 
For 

m = p 3 \t + /? 4 r +1/2 cos r, /? 2 sin r, ft ) T 

wrai/i r = xln |i + 04\ + (5 and /3 a , f3^ = const, where 0\ + /? 2 = 1, i^e operator D + 2/?4<9 t + 2xrJi 2 
induces the operator 

D\ + 2/3 4 <9 y3 - 2/3 1 xJ 1 1 2 + 2av 3 d v3 , 

where D\ = yid Vi + 2y 3 d y3 — v % d v i — 2qd q , if the vector-functions h l are chosen in form $Z.1\ ). In 
all other cases the basis elements of the MIA of (2.S) are not induced by operators from A(NS). 



Note 2.8 The invariance algebras of systems of form (2J>) with different parameter-functions p 3 = 
p 3 (t) and p 3 = p 3 (t) are similar . It suggests that there exists a local transformation of variables 
which make p 3 vanish. So, let us transform variables in the following way: 

v l = (V + i yi/ 9 3 (t))e-^W, v 3 = v 3 , (2.12) 
q = qe -P(t) + § wl/i ((p3( f )2) _ 2 p 3 (t))e-^\ 
As a result, we obtain the system 

v\ + &v % j - Vjj + qi + p l {yz)v^ = 0, 
vl + v j v 3 - v 3 j = 0, 
5f = 

for the functions v a = v a (yi, y 2 , 2/3) and q = (/(yi, 2/2) 2/3)- Here subscripts 1, 2, and 3 denote 
differentiation with respect to y\, y 2 , and 2/3, accordingly. Also p % (y 3 ) = p t (t)e~2 p ^ . 
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3 Reduction of the Navier-Stokes equations to systems of PDEs 
in two independent variables 



3.1 Ansatzes of codimension two 

In this subsection we give ansatzes that reduce the NSEs to systems of PDEs in two independent 



variables. The ansatzes are constructed with the subalgebrical analysis of A(NS) ( see Subsec. [A.3 
) by means of the method discribed in Sec. |B] . 



1. u 1 = (rR) 1 ((x\ — kx2)w 1 — X2W 2 + x\x^r l w 3 ), 
(rR)~ 1 ((x2 + xxijw 1 + xiw 2 + X2X3r~ 1 w 3 ), 



V 2 



u 3 = x 3 (rR) V - R 1 w 3 , 
V = R~ 2 s, 



where Z\ = arctana^/^i — xlnR, Z2 = arctanr/x3, k > 0. 

Here and below w a = w a (zi, Z2), s = s(z±, z 2 ), r = (x 2 + x 2 ) 1 ^ 2 , R = {x\ + x\ + x 2 ) 1 / 2 , 
x, e, a, n, and v are real constants. 



2. u 1 = \t\ 1 / 2 r 1 (xiw 1 — X2W 2 ) + l x\+x\r 2 , 

-2 

(3.2) 

|i| 1 s — \r~ 2 + \t~ 2 R 2 + e\t\~ x arctanx 2 /xi, 



V? = \t\ X I 2 T 1 (X2W 1 + X\W 2 ) + \t 1 X2+X2T 



U 3 = \t\ 1 / 2 W 3 + XT l W 2 + \t 1 X%, 



where z\ = \t\ 1 / 2 r, z 2 = \t\ 1 ' 2 xs — ^arctan X2/X1, x>0, e > 0. 

3. u 1 = r^ixiw 1 — X2W 2 ) + xir~ 2 , 
u 2 = r~ 1 (x2W 1 + x\w 2 ) + X2r~ 2 , 
u 3 = w 3 + xr~ 1 w 2 , 
p = s — \r~ 2 + e arctan x 2 /xi, 
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(3.3) 



where z\ = r, Z2 = X3 — xarctanx2/xi, x £ {0; 1}, £ > if x = 1 and e G {0; 1} if 
x = 0. 

4. u 1 = |t| _1 / 2 (^ti; 1 + zvti; 3 ) cost - Itj" 1 / 2 ^ 2 sin r+ 

+z^t _1 cost + ^t _1 xi — >rf -1 X2, 
w 2 = |t| _1 / 2 (^ 1 + uw 3 ) sinr + 1 1| ~ 1 / 2 w 2 cosr+ 

+z/£t~ 1 sin t + \t~ x x 2 + xt _1 xi, (3.4) 
u 3 = |t|- 1 / 2 (— i^w 1 + uw 3 ) + //£t _1 + 5t _1 x 3 , 
p = |t|-i a - lf 2 e + |r 2 i? 2 + ^t- 2 r 2 + 

+e\t\~ 3 l 2 {vxi cost + sinr + fix^), 

where 

= |t| _1//2 (/xxi cost + ^2 sinr — 1/X3), 
z 2 = |t|~ 1 ^ 2 (x2 cost — xi sinr), 

£ = a{yx\ cost + 2^X2 sinr + /XX3) + 2kv(x 2 cost — xi sinr), 
T = x\n\t\, k>0, n>0, v > 0, // 2 + i/ 2 = l, cre = 0, e > 0. 

5. u 1 = {tl-^w 1 + ^ _1 xi, 

n 2 = |t|-V2 w 2 + l t -l X2) 

p = |t|-i a _ 1^-23.2 + |r 2 i? 2 + £ |t|- 3 / 2 x 3 , 

where 

^1 = |i|~ 1/2 aTi, z 2 = |tr 1/2 x 2 , cre = 0, e>0. 

6. -u 1 = (/Uti; 1 + uw 3 ) cost — w 2 sin i + z/£ cost — X2, 



(3.5) 



« 2 



(/Ltti; 1 + uw 3 ) sint + w 2 cost + z/£sint + xi, 
-u 3 = (— uw 1 + fiW 3 ) + (J,£, 

p = s — 7j£ 2 + ^r 2 + e(yx\ cost + z^X2 sin t + /UX 3 ), 

where 

zi = {pLX\ cos t + [1x2 sint — 2^x3), 
^2 = (X2 cos t — xi sint), 

£ = cr(^xi cost + z/X2 sint + /iX3) + 2i/(x2 cost — xi sint), 
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(3.6) 



fi > 0, v > 0, ji/ 2 + z/ 2 = 1, ae = 0, e > 0. 



7. -u 1 = u; 1 , u 2 = w 2 , u 3 = u> 3 + axs, 



p = S — t;(J 2 X3 + EX3, 



where 



zi = xi, z 2 = x 2 , as = 0, eG{0;l}. 

5. U 1 = XlW 1 — X2T~ 2 (w 2 — x(t)), 

u 2 = x 2 w l + x±r~ 2 (w 2 — 

"W 3 = (/f(i)) _1 (w 3 + Pt(t)x3 + earctanx2/xi), 

p = s - 5/9 tt (t)(/9(t)) _1 x| + xt(i)arctanx 2 /xi, 



where 



Zl =t, z 2 = r, £G{0;1}, X , P € C°°((to, 



9. u = w + A 1 (n*-x)m^-A 1 {k-x)k t , 

p = s — \\~ l (m l tt ■ x)(n l ■ x) - \\~ 2 {m\ t ■ k){n % ■ x)(k ■ x), 



where 



zi=t, z 2 = (k-x), m l e C°°((to,ti),R 3 ), 

I n -| n — '101 O — ' 

fh tt ■ m — m ■ fh tt = 0, k = in x m , n = m x k, 
n 2 = kxm 1 , X = \(t) = k-k^0 Vt€(t ,*i). 



(3.7) 



(3.8) 



(3.9) 



12 



3.2 Reduced systems 



Substituting ansatzes (3.1)— (3.9) into the NSEs (1.1), we obtain the following systems of reduced 
equations: 

1. W 2 w\ + W 3 W 2 — W lr W 3 COt Z 2 — (w 1 ) 2 — (iV 2 + KW 1 ) 2 sin 2 Z 2 — 



-(w 3 ) 2 - {{k 2 + sin 2 z 2 )w\ 1 + w\ 2 



2w\ 



2w 2 - 



—2w l ) sin Z2 + w 2 cos z 2 — w 1 sin 1 z 2 — (2s + xs\) sin 2 z 2 = 0, 

w 2 w 2 + w 3 w 2 + w 3 (w 2 + 2kw v ) cot z 2 — 

-^((w 1 ) 2 + (w 3 ) 2 + (w 2 + yew 1 ) 2 sin 2 z 2 )- 

-((k 2 + sin~ 2 z^w^ + w 2 2 + 3>m; 2 + 2>fr(u;2 + xw\ + w 1 ))- 



■ sin z 2 + (2w;{ + 2iuf cot z 2 



w 



2>rtt; 1 )sin 1 z 2 - 



-(w 2 + 2xw 2 ) cos z 2 + 2><s sin 2 Z2 + (1 + x 2 sin 2 z 2 )s\ = 0, 



w 2 w 3 + w 3 w 3 



w 3 ) 2 cot Z2 — {w 2 + ot 1 ) 2 sin 22 COSZ2- 
— ((x 2 + sin -2 ^2)^11 + ^22 + ^1 + 2^2) sinz2 + 
+(2 , u; 1 + uif + w\ + xitfj) cos z 2 + S2 sin 2 Z2 = 0, 



w 1 + w 2 + w' 2 



0. 



(3.10) 



Hereafter numeration of the reduced systems corresponds to that of the ansatzes in Subsec. 
Subscripts 1 and 2 denote differentiation with respect to the variables z\ and z 2 , accordingly. 



2-3. w l w\ + w 3 w\ - z x 1 w 2 w 2 - [w\-y + (1 + * 2 z 1 2 )w 22 )- 
-2xz x 2 w 2 + si = 0, 

w l w 2 + w 3 w 2 + z^w^w 2 — (w 2 ! + (1 + k 2 z 1 ~ 2 )w 22 ) + 
+2xz 1 ~ 2 w\ + 2z 1 ~ 2 w 2 - xz l T l s 2 + ez{ 1 = 0, 

w 1 w 3 + w 3 w 3 — 2xz l ~ 2 w 1 w 2 — (wfi + (1 + m 2 z x 2 )w 22 )+ 
+2k(z^ 2 w 2 )i - 2>?z{ 3 w\ + (1 + x 2 z^ 2 )s 2 - exz{ 2 = 0, 

w\ + w 2 + z^w 1 +7 = 0, 



(3.11) 



where 7 = ±3/2 for ansatz ( |3.2[ ) and 7 = for ansatz (3.3). Here and below the upper and lower 
sign in the symbols "±" and "=p" are associated with t > and t < 0, respectively. 
4-7. For ansatzes (|3.4j) — (|3TT[) the reduced equations can be written in the form 



w wt 



w}j + s\ + a 2 w 2 



0. 



w l w 2 — w\ + s 2 — a 2 w l + a\w 3 = 0, 
w l w 3 — wft + a^w 3 + 05 = 0, 



(3.12) 



w\ 



a 3 
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where the constants a n (n = 1, 5), take on the values 



4. a.\ = ±2xi/, «2 = T2x/x, 03 = =f(<t + 3/2), 04 = ±<r, Q5 = e. 

5. ai = 0, a2 = 0, 03 = =f(o" + 3/2), 04 = ±cr, 05 = 6. 

6. ai = 2z/, 02 = — 2/U, 03 = — <7, 04 = <J, 05 = £■ 

7. ol\ = 0, «2 = 0, «3 = — it, 04 = <J, 05 = e. 

8. V)\ + (w 1 ) 2 — -Z 2 " 4 ( u ' 2 ~~ X) 2 + Z2W l w\ — W22~ 

-3z 2 Wr, + z 2 ~ 1 s 2 = 0, 

o 19 9 19 

tl^ + Z2W ^2 ~~ "^22 + z 2 ^2 = 0> 

+ ^W 1 ^ - U>22 - Z 2~ 1,W 2 + ^2~ 2 (^ 2 - X) = 0, 
2U) 1 + Z 2 ^2 + pi/p = 0. 

9. uJi — Au?22 + s 2 k + A _1 (rf • u;)m| + ^2? = 0, 
k ■ W2 = 0, 
where y\ = t and 

e = e(i) = 2A~ 2 (m f 1 • rh 2 - rh 1 ■ fh 2 )k t x fe + X~ 2 {2k t ■ k t - ht • 
Let us study symmetry properties of reduced systems fl3.10p and fl3.ll ) 

Theorem 3.1 T/ie ML4 of \3.1C ) is given by the algebra < d\ >. 



(3.13) 

(3.14) 
(3.15) 
(3.16) 
(3.17) 
(3.18) 



Theorem 3.2 The MIA of (3.11) is given by the following algebras: 

a) < 82, d s ,D 2 = Zidi — w a d w a — 2sd s > if 7 = x = e = 0; 

b) <d 2 ,d s > if ( 7 ,x,e) 7^(0,0,0). 



All the Lie symmetry operators of systems ( |3.1C ) and ( |3.11| ) are induced by elements of A(NS). 
So, for system ( |3.10| ) the operator d\ is induced by J\2- For system ( |3.11| ), when 7 = (7 = ±3/2), 
the operators D 2 , 82, and d s (82 and d s ) are induced by D, R(0, 0, 1), and Z(l) (i?(0, 0, |t| -1 / 2 ) 
and accordingly. Therefore, the Lie reductions of systems ( |3.1C| ) and ( |3.11| ) give only 

solutions that can be obtained by reducing the NSEs with three-dimensional subalgebras of A(NS) 
immediately to ODEs. 

Investigation of reduced systems (fU3|)-(fl|), and (|Q3) 

is given in Sec. [| and 



4 Reduction of the Navier-Stokes equations to ordinary differen- 
tial equations 
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4.1 Ansatzes of codimension three 

By means of subalgebraic analysis of A(NS) (see Subsec. |A.3|) and the method described in Sec. || 
one can obtain the following ansatzes that reduce the NSEs to ODEs: 

1. u 1 = x\R~ 2 (p l — X2(Rr)~ 1 f 2 + xix%r~ l R~ 2 tp 3 , 

u 2 = X2R~ 2 (f 1 + xi(Rr)~ 1 (p 2 + X2Xsr~ 1 R~ 2 (p 3 , 

u 3 = x^R~ 2 ip l — rR~ 2 ip 3 , 

p = R- 2 h, 



u 2 = X2<£ 1 + X\T 2 ip 2 , 



u = (f + crx;j + v arctan x-^/xi, 

p = h — ^a 2 x 2 + ei arctan x^fxx + £2^3, 
where u = r, va = 0, £20" = 0, and for a = one of the conditions 

v = 1, £1 > 0; v = 0, £1 = 1, e 2 > 0; v = £1 = 0, £ 2 G {0; 1} 
is satisfied. 

Two ansatzes are described better in the following way: 
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(4.1) 



(4.2) 



where u = arctan r/x^. Here and below (p a = (p a (oj), h = h(uj), 
r = (x\ + xl) 1 / 2 , R = (xl + xl + xl) 1 / 2 . 

2. u 1 = r^ 2 (x\ip 1 — X2(f 2 ), u 2 = r~ 2 (x2(p 1 + x\{p 2 ), 
u 3 = r _1 (^ 3 , p = r~ 2 h, 

where uj = arctanx2/xi — >rlnr, >c > 0. 

3. u 1 = xi\t\~ l ip l - X2r~ 2 tp 2 + \x\t~ l , 
u 2 = x 2 \t\- l ip 1 + xir~ 2 tp 2 + ±x 2 t~ 1 , 

u 3 = \t\-V 2 ( p3 + (ff+i)x 3 t- 1 + v\t\ 1 /H- 1 ax<Afmx 2 /x 1 , (4.3) 
p = \ t \-^h+\t- 2 R 2 -\a 2 xlt- 2 + 

+£i|t| _1 arctanx2/xi + ^2^3 1 ^| 3 / 2 , 
where uj = \t\~ 1 / 2 r, va = 0, £2<r = 0, e\ > 0, v > 0. 

4. it 1 = x\tp l — X2r~ 2 (p 2 , 



(4.4) 



5. The expressions for u a and p are determined by (|2.l| ), where 



v 1 = anp 1 + a 2 f 3 + buUi, 

v 2 = Lp 2 + b 2 iUJi, 

v 3 = a 2 ip x - anp 3 + b 3i Ui, 

P = h + CuUli + C2iU!U>i + \dijOJibJj. 



(4.5) 



(4.6) 



In formulas (4.5) we use the following definitions: 

^i = aiyx + a 2 y 3 , u 2 = y 2 , w = u 3 = a 2 y\ - a\y 3 ; 
ai = const, a\ + a 2 = 1; a 2 = if 71 = 0; 

7! = -2x, 72 = — § if i > and 71 = 2x, 72 = | if i < 0. 
6„i, fij, Cy, and dij are real constants that satisfy the equations 
hi = aiBi, b 3i = a 2 Bi, c 2i + a 2 7i&2i = 0, 

621 Bi + b 22 b 2i - 7iai-Bj + d 2i = 0, 
BxBi + £ 2 & 2 j + 7iai-Bi + d u = 0, 

(Si + b 22 )(B 2 + ai7i - 621) = 0. 

6. The expressions for u a and p have form Q2.2j ), where u a and 5 are determined by ( |4.5| ), (|4.G|). 

and 71 = — 2x, 72 = 0. 

Note 4.1 Formulas faj.3i) and (4-t) determine an ansatz for system (2.7), where equations 

are the necessary and sufficient condition to reduce system (2.1) by means of an ansatz of form 

7. u 1 = ip 1 COSX3/?? 3 — if 2 sm.x 3 /rf' + x\0 (t) + x 2 8 2 {t), 
u 2 = ip l sina^/^ 3 + ip 2 cosx 3 /rj 3 — x\6 2 {t) + x 2 6 l (t), 
u 3 = ip 3 + ri 3 (rf i )~ 1 x 3 , 

P = h - IvUv 3 )" 1 ^ ~ Wtt r l j {v i ri i )~ 1 r 2 , 
where u> = t, 



(4.7) 



r] a G C°°((to,*i),K), r/VO, rfrf^Q, rtW - rfrj 2 G {0; ±}, 



0i = 4rf{rfr?r\ 2 = (njn 2 - ^Hr/ty')- 1 . 



P 



h — A (m" t • x)(n a • x)+ 



(4.8) 



-fiA 2 {m\ t ■ m a ){n a ■ x){n h ■ x), 
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where u = t, rn a G C°°((t , h), E), m£ ■ fh b - fh a ■ m\ t = 0, 
A = A(t) = (m 1 x m 2 ) • m 3 ^ Vie(io,*i), 

->1 ->2 -» 3 -^2 -> 3 -> 1 ->3 -» 1 ^2 

n = m x m , n = m x m , n = m x m . 
4.2 Reduced systems 

Substituting the ansatzes 1-8 into the NSEs fll.l| ), we obtain the following systems of ODE in the 
functions ip a and h: 

1. </>Vi - <^V a - y>L - </>i cot u - 2h = o, 

+ ¥> V cot w - yl w - ipl cot uj + ip 2 sin" 2 w = 0, 

tp 3 ip 3 - ip 2 ip 2 cot w - tp^ - tpZ cot lu + ip 3 sin -2 w+ ( 4-9 ) 
-2p* + hu = 0, 

V 1 + V? 3 + ^ 3 cot w = 0. 

2. (c^ 2 - V)¥>i - (1 + x 2 )^ - V - <^V - x^, -2/i = 0, 
(y> 2 - V)</£ - (1 + - 2(^ 2 + p£ ) + K = 0, 

(4.10) 

((^ 2 - x^ 1 )^ _ (1 + ^)<^ _ ^3 _ ^3 _ 2 ^3 = Qj 

vl - *pI = o. 

3-4. <^v - w"VV 2 + w^Vi - - 3w-Vi + = o, 
<^VS - + w ~V 2 + ei = o, 

(4.11) 

uip l ifl + criv? 3 + vu> V -<plu-u) V 3 +^2 = 0, 
lip 1 + + a 2 = 0, 

where 

3. ai = a, a 2 = (<x + §) if t > 0, 
cji = -a, a 2 = -(a + §) if t < 0. 

4. O"! = 0"2 = o. 

5-6. (/^Vi, - ¥>L - Mii</> 1 + en + c 2 iw = 0, 

- <#L ~ ^2HP 1 + C12 + C 2 2W + 72a 2 V? 3 = 0, 

(4.12) 

- <fitu + 71 <W 2 + /la; = 0, 

ft = <r, 
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where /in = —B\, /ti 2 



7i a i> M21 



7. + #V + v - (v 3 )- W + fa 3 rV 



&21 + 7l a l, M22 

0, 



-622, <r = 71 - J3i 



->22- 



95; 



2,„2 



0, 



^ + 77f(r? 3 ) V 3 = 0, 



(4.13) 



0. 



29 1 +i]f(7] 3 )- 1 = 0. 
8. u) + \- l {n b -p)m b t 
n a ■ mf = 0. 

4.3 Exact solutions of the reduced systems 



(4.14) 



1. Ansatz (4.1) and system ( |4.9| ) determine the class of solutions of the NSEs (1.1) that are called 
the steady axially symmetric conically similar flows of a viscous fluid in hydrodynamics. This class 
of solutions was studied in a number of works (for example, see references in ||lq| ). For ip 2 = it was 
shown, by N.A.Slezkin |54|, that system Q4.9| ) is reduced to a Riccati equation. The general solution 
of this equation was expressed in terms of hyper geometric functions. Later similar calculations were 
made by V.I.Yatseev [^] and H.B. Squire |35|. The particular case in the class of solutions with 



■ is formed by the Landau jets |24J. For swirling flows, where ip 7^ 0, the order of system 
can be reduced too. For example 



(f 2 (p 2 sin 2 u) — sinu;^^ sin 1 u) w + 2$^ cot uj + 2<3? 



an arbitrary solution of ( |4.9| ) satisfies the equation 
const, 



(tp^ — ^ip 3 tp 3 ) sin 2 uj — ip 3 coscjsinw, and the Yatseev results [38] are completely extended 



where 3> 

to the case ip 2 sinco = const. 
2. System (4.10) implies that 



<p 2 = xf 1 + Ci, 

h = x{l+ k 2 )^ + {2^ + 2- xC^ip 1 + C 2 , 

(1 + x 2 )^ + (4x - d)(pi + 9? V + V+ 
+ (l + ^ 2 )- 1 (C 1 2 + 2C 2 ) = 0, 



(4.15) 



(1 + - 2 ) 



(Ci-2x)<fi, + {l + <p 1 ) v ? = 0. 



If ip 3 = 0, the solution determined by ansatz ( |4. 10|) and formulas (4.15) coincides with the Hamel 
solution |l8|, |23|. In Sec. ^ we consider system (|6.14 ) which is more general than system (4.1C). 
3-4. Let us integrate the last equation of system ( [4.11 ), i.e., 



2°2- 



(4.16) 



Taking into account the integration result, the other equations of system (4.11) can be written in 
the form 



hu 



1 2 



((Ci + l^- 1 - \a 2 ^l 



£1, 
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<fL> - ((^i - x ) w 1 - - a ^ = vuj V + £2. 



(4.17) 



Therefore, 

h = J u!~ 3 (p 2 (p 2 dio - \C 2 uj~ 2 - gcrfo; 2 , 

ip 2 = C 2 + C 3 J \Lo\ c ^ +1 e-^ 2 duo+ 

+ei / M^ +1 e-3 ff2aj2 (/ Iwl-^-^i^ 2 ^)^. 
If o"i = 0, it follows that 

<p 3 = C 4 + C 5 J \uj\ c ^- 1 e-^ 2 duj+ 

+ j\uf^- 1 e-^ 2 (j \to\- Cl+1 e^ u2 {e 2 + uu- 2 <f?)dw)du. 



(4.18) 
(4.19) 



(4.20) 



Let di / (and, therefore, v = 0). Then, if 02 7^ 0, the general solution of equation ( 4.171 ) is 
expressed in terms of Whittaker functions: 

^3 = \^-l e -\^ w{ _ aia -l + 1^ _ l f l Cu 1^2^ 

where W{x, fj,, r) is the general solution of the Whittaker equation 

4r 2 W TT = (r 2 - Axt + Afi 2 - l)W. (4.21) 
If a 2 = 0, the general solution of equation ( 4.16| ) is expressed in terms of Bessel functions: 



w|3 c iZi c ((-ai) 1 /^), 



V2, 



where Z v {r) is the general solution of the Bessel equation 
t 2 Z tt + rZ r + (r 2 - z^ 2 )Z = 0. 



(4.22) 



Note 4.2 // o"2 = 0, all quadratures in formulas M-.18Q -(4-2C ) are easily integrated. For example 



C 2 + C 3 ln\uj\ + \eioJ 2 
C 2 + C 3 \uj 2 + ^eiw 2 (lnw - \) 
C 2 + C 3 (C 1 + 2)~ 1 \uj\ c ^ +2 -\e x C^uj 2 if Ci^-2,0. 



Ci = -2, 
d = 0, 



5-6. Let <7 = 0. Then the last equation of system ( 4.12j ) implies that c^ 3 = Co = const. The 
other equations of system fl4.12j ) can be written in the form 



h = —Jia 2 J ip 2 (uj)du!, 



(4.23) 



where v\\ = en, v 2 \ = c 2 \, v\ 2 = c^+^^Cq, v 22 = c 22 . System ( |4.23j ) is a linear nonhomogeneous 
system of ODEs with constant coefficients. The form of its general solution depends on the Jordan 
form of the matrix M = {fiij}. Now let us transform the dependent variables 

if 1 = eijijj j , 
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where the constants e« are determined by means of the system of linear algebraic equations 

^ijfijk — t^ij&jk (£) ji k — 1, 2) 

with the condition det{ejj} 7^ 0. Here M = {/%} is the real Jordan form of the matrix M. The 
new unknown functions ip % have to satisfy the following system 

V>L - CoC + = vu + "2iU, (4.24) 

where = eij&ij, u 2 i = &ijV 2 j- Depending on the form of M, we consider the following cases: 
A. detM = (this is equivalent to the condition detM = ). 

i. M = ( jj e Q ^j , where e G {0; 1}. Then 

^ = d + C 2 e c ^ - ^ 22 C W - (P 12 - ^Cq- 1 )^ V (4.25) 
V' 1 = C 3 + C 4 e c ^ - \v 21 C^u? - {u u - 9 21 C 1 )C 1 uj+ 

+e(-^ 22 C -V - \(v X2 - 29 22 C 1 )C 2 u; 2 + 

+(Ci + (£ 2 i - 2P 22 C - 1 )C ( 7 2 )C ( 7 1 ^ - C 2 Co 1 Lue c ° u 
for Co 7^ 0, and 

^,2 = Cl + C 2 u; + ±£ 22 u; 3 + \ vi 2 u 2 , (4.26) 
^ = C 3 + C 4 w + i(£*n - C 2 )w 3 + ±(% - d)u> 2 - 



for C = 0. 



ii. M = ( ^ y ' w ^ ere Hx ^ ^\{0}- Then the form of ip 2 is given either by formula 
(H21) for C / or by formula (gjjf ) for C = 0. The form of ^ is given by formula (|4.28|) (see 



below) . 

B. detM / (this is equivalent to the condition detM 7^ 0). 

i. M = ( ° ) , where *g G M\{0}. Then 

V u ^ 2 y 

?/, 2 = ^^-^ + (z>i 2 - CqD 22 ^ 1 )^ + Ci# 21 (w) + C 2 22 (uj), (4.27) 

= ^ix^o + (hi - Co^i^ 1 )^ 1 + C 3 8 n (u) + C 4 # 12 (c;), (4.28) 

where 

20 



O il (u) = exp(i(C - VA)cu), 9 i2 {uo) = exp(±(C + VA)w) 
if Di = C 2 - Axi > 0, 

9 a (u) = e5 Cow cos(±v/=Aw), i2 {uj) = e^ Co " sm(±^D~iLj) 
if Di < 0, 

0*1(0;) = e^ CoU , e i2 (uj) = coe^ CouJ 
if Di = 0. 



ii. M = [ X2 ^ j , where X2 G R\{0}. Then the form of i)j 2 is given by formula ( 4.27| ), 



X2 

and 

Ip 1 = (pll - (p\i - C i>22^ 1 ) x 2 1 ~~ ^0(^21 - ) x 2 ) x 2 X + 

+(P 21 - z> 22 ^ 1 )x 2 x cj + C 3 fl 21 (w) + C 4 fl 22 (u;) - CiT/*^), 

where 

^'(u;) = J D 2 - 1 u;(202i(u;)-Co0^(a;)) if £> 2 / 0, 
r ] 1 {uj) = lu 2 e^ Cou; , r?{Lo) = \uj^e^ Cou) if D 2 = 0. 

iii. M = [ Xl I , where xi G R, x 2 / 0. Then 

^)1 = {xi»li)~ X {v%\X\ + Z>22^2)w + (xiXiY 1 (Vuxi + &12X2)- 

-C (xiXi)- 2 (i>2l{x2 ~ *?) - i>222^lX 2 ) + C n ln (w), 
V> 2 = (^^) _1 (-Z>21X!2 + i>22^l)w + (xi^) _1 (-Z>n^2 + ^12^1)- 

-C (^^) _2 (z : '2i2xi^2 + ^22(^2 - *l)) + C n 9 2n (uj), 
where n = 1, 4, 



7 = V'(C 2 -4>, 1 )2 + (4x2) 2 , 



A = £^2(7 + C 2 -4*0, /3 2 = ^^2(7- C 2 + 4*0, 

= 22 ( w ) = exp((iCo-#Ow)cos/92W, 
-0 21 H = 12 (w) = exp((iC -fr)w)sin#jw, 
0i3( w ) = ^24^) = exp((iC + pju) cos 
023( w ) = _ i4( w j = exp((iCo + AV)sm/32^- 
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If a 0, the last equation of system ( 4.12 ) implies that ip 3 = ato (translating u>, the integration 



constant can be made to vanish). The other equations of system ( [hl2|) can be written in the form 



h = — 7iQ2 / <p 2 (u})du — \a 2 ui 2 , 



flu ~ <*u(pl + fiijtp 1 = v u + u 2i uJ, (4.29) 

where v\\ = c\\, v 2 \ = c 2 i, v\ 2 = c\ 2 , v 22 = c 22 + 72020". The form of the general solution of 
system ( 4.29| ) depends on the Jordan form of the matrix M = {fiij}. Now, let us transform the 



dependent variables 

ip l = eijip j , 

where the constants are determined by means of the system of linear algebraic equations 

&ijfl'jk — l^ijGjk (J'tji k — 1, 2) 

with the condition detje^'} 7^ 0. Here M = {/%} is the real Jordan form of the matrix M. The 
new unknown functions ip l have to satisfy the following system 

i>lw - vu^l, + ilijip j = vu + D 2i io, (4.30) 
where uu = eij&ij, v 2 i = £ijV 2 j- Depending on the form of M, we consider the following cases: 
A. detM = (this is equivalent to the condition detM = ). 

i. M = ( jj £ Q ^j , where e e {0; 1}. Then 

i/j 2 = C 1 + C 2 J ei aw ' 2 duj - a- 1 9 22 u+ 

(4.31) 

+D l2 / e^ 2 (J e-^ 2 du)duj, 
^1 = C 3 + C 4 J eh™ 2 6w - a- 1 9 21 oj+ 

+ Je^ auj2 (Je-^ 2 (i)n-eiP 2 )duj)dio. 



ii. M = ( a ] . Then the form of tp 2 is given by formula ( 4.31| ), and 



v / 

V/ = C 3 lu + C\(uj f e^dio - o- l eh™ 2 ) + ff- 1 v n + 
+a- 1 D 21 {auj J e^ 2 X 1 (u)dw - e^™ 2 A 1 (u>)) , 
where A x (w) = J e'^div. 



hi. M 



^ ^ ' wnere x i e K\{0;cj}. Then ip 2 is determited by ( [4.31 ), and the form 
(P~3l ) (see below). 



of V' 1 is given by fl4.33| ) (see below) 
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B. det M^O, det{flij—a5ij} = (this is equivalent to the conditions det M^O, detj/Xjj — <x<%} = 
0; here 5ij is the Kronecker symbol). 

i. M = ^ g ^ j , where e G {0; 1}. Then 

V; 2 = Ciw + C 2 (w J e^duj - ff^eW) + (j- 1 z>i2 + 

(4.32) 

+ C j- 1 P 22 (cJw/e5 <J - 2 A 1 H^ - e^A^w)), 
^1 = C 3 u; + C 4 (^ J e^ 2 dw - o~ x e^ 2 ) + a^v^ 

+aujfe 1 2™ 2 \ 2 (u;)doj - e^ 2 \ 2 {uj)+a- l {D 2 iuj - e^ 2 ), 
where A 1 (w) = Je"^ 2 ^, A 2 (w) = a' 1 f e~l au3 \v 2 i - eV£)dw. 

ii. M = ^ ^ ^ J , where xi G K\{0;<r}. In this case V> 2 is determined by ( 4.32|) , and the 
form of tp 1 is given by ( |4.33| ) (see below). 

C. detM 7^ 0, det{flij — a5ij} ^ (this is equivalent to the condition detM / 0, det{fiij — adij} ^ 
0: here 5ij is the Kronecker symbol). 

i. M = [ *J ° ) , where xg G R\{0;cr}. Then 

iP 1 = x^vn + (x! - o-) _1 ^2i^ + M _1 / 2 e3 CTa;2 - 
■(CaM^xiff- 1 + ±, \, \au?) + C 4 M(ix lC 7- 1 + \, 
= >q x v X i + (x 2 - o-)- 1 ^^ + |w|- 1 / 2 ei CT ^ 2 - 

(4.34) 

•(CiM(ix 2 a- 1 + ±, ±, \au?) + CW^a" 1 + |, -|, ±^ 2 )), 

where M(x, fj,, r) is the Whittaker function: 

M(x,/x,r) =r^ + ^e-5 T iF 1 (i+ / u-x,2/i+l,r), (4.35) 

and i-Fi(a, 6, r) is the degenerate hypergeometric function defined by means of the series: 

^ q(a + l)...(a + ra- 1) r n 
lFl(o ' 6 ' T) = 1 + g 1 6(6+l)...(6 + n-l) 



(4.33) 

1,1 11. 



6/0,-1,-2, 



ii. M = ( Xl ^ 2 , where *g el, x 2 ^ 0. Then 
y x 2 xi y 

^ x = (xjXj) _1 (xiz>n + x 2 z>l 2 ) + 

+ ((xi - a) 2 + x 2 .) -1 ^! - <r)i>2i + x 2 i> 22 )u;+ 
+C 1 Rer ] 1 {uj) - C 2 Imrj 1 (to) + C 3 Rer] 2 (uj) - C 4 Imrj 2 (to), 
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if 2 = (x i x i )" 1 (-x 2 z>n + >tj>vz)+ 

- a) 2 + x|) _1 (— x 2 z>2i + (x\ - cr)h>22)u+ 
+C 1 Imi ] 1 {uj) + C 2 Rerj l (Lo) + Cslmr? 2 ^) + C 4 Rerj 2 (Lo), 

where 

r? 1 ^) = M(|(^ + x 2 \)a- 1 + ±, |, \<jlo 2 ), 

r, 2 {u) = M(±(*i + ^cj- 1 + ±, -±, ±c^ 2 ), i 2 = -1. 



iii. M = ( ^ ) , where X2 £ R\{0;er}. Here the form of ip 2 is given by (133), and 

V ;1 = (^11 - £l2><2~ ) x 2 + (^21 - ^22^2 - ^)~ 1 )(X2 - ay 1 UJ+ 

^uj\~ l ' 2 e^ 2 (Cz6\t) + C a 6 2 {t) - g- 1 6\t) $ T~ l 6 2 {T)C i e i (r)dT+ 



-a 



where r = ^auj 2 , 

6 l {r) = Md^a- 1 + i, i,r), fl 2 (r) = M^^- 1 + \, -\,r). 

Note 4.3 TTie general solution of the equation 

ipuu - ffw^ - (n + l)o-ip = 0, 
where n is an integer and n > 0, is determined by the formula 
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Note 4.4 If function tj) satisfies the equation 

ipajui — (TWipw + = (x^— cr), 
i/ien / ^(uj)duj = (x+ c) -1 (<7u;?/> ~~ + Ci- 



7. The last equation of system (4.13) is the compatibility condition of the NSEs ( |1.1| ) and ansatz 



L7| ). Integrating this equation, we obtain that 
r, 3 = CoW)' 1 , Co 7^0. 
As = — f? 3 (f? 3 ) _: V 3 = 20 1 ( / 9 3 , y? 3 = C^rfn 1 . Then system (433) is reduced to the equations 

<pi = xW -x 2 (^V 2 , 

(4.36) 

v? 2 = xW 1 + xW 2 , 
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where = -C^irfrf) 2 - 6 1 and X 2 = 2 ~ C^C^^rf) 2 . System (gjg ) implies that 
Lp 1 = exp(/x 1 (^)^)fCicos(/x 2 (^)^) - C 2 sin(/x 2 (^)^)), 



V? 2 = exp(fx 1 (u)dv)(c 1 sm(fx 2 (u)duj) + C 2 cos(J x 2 (^)dto] 



8. Let us apply the trasformation generated by the operator R(k(t)), where 
kt = \- 1 (n b ■ k)m\ - <p, 



to ansatz ( |4.8[) . As a result we obtain an ansatz of the same form, where the functions (p and h are 
replaced by the new functions <p and h: 

(p = <p — A~ 1 (n a • fc)m£ + fc t = 0, 

h = h- \~ l (fh u • fc)(n a • jfc) + ^A- 2 (m^ • m a )(n a • k)(n b ■ k). 

Let us make h vanish by means of the transformation generated by the operator Z(—h(t)). There- 
fore, the functions <pf and h can be considered to vanish. The equation (n a • mf) = is the 



compatibility condition of ansatz (4.8) and the NSEs ( |1 . 1|) . 

Note 4.5 The solutions of the NSEs obtained by means of ansatzes 5-8 are equivalent to either 



solutions (5.1) or solutions (5.5). 



5 Reduction of the Navier-Stokes equations to linear systems of 
PDEs 



Let us show that non-linear systems 8 and 9, from Subsec. 3.2, are reduced to linear systems of 
PDEs. 

5.1 Investigation of system (|3Tl7|) - (|3TT8| ) 



Consider system 9 from Subsec. [T^, i.e., equations J3.17 ) and ( 3.1 8| ) . Equation ( 3.18 ) integrates 
with respect to z 2 to the following expression: 

k ■ w = ip(t). 

Here i/j = ip(t) is an arbitrary smooth function of Z\ = t. Let us make the transformation from the 
symmetry group of the NSEs: 

u(t, x) = u(t, x — l) + lt(t), 

p(t, x) = p(t, x-l) - l tt (t) ■ x, 

where la • m ! — I ■ rh\ t = and 

k-(h- A _1 (n i • l)mi + A~ x (fc • l)k t ) + ip = 0. 



This transformation does not modify ansatz ( |3.9| ), but it makes the function tp(t) vanish, i.e., 
k ■ w = 0. Therefore, without loss of generality we may assume, at once, that k ■ w = 0. 
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Let f 1 = f l (zi, z 2 ) = m 1 -w. Since m\ t ■ fh 2 — fh 1 ■ fh 2 t = 0, it follows that rh\ ■ fh 2 — fh 1 ■ fh 2 = 
C = const. Let us multiply the scalar equation (3.17) by m 1 and k. As a result we obtain the linear 
system of PDEs with variable coefficients in the functions f 1 and s: 

fi ~ x fL + CX^dm 1 ■ fh 2 )/ 1 - (ttV • m 1 )/ 2 )- 

-2C\- 2 ((k x k t ) ■ m l > 2 = 0, 

s 2 = 2A- 2 K • k t )f + \-\ht -k-2k t - k t )z 2 . 

Consider two possible cases. 

A. Let C = 0. Then there exist functions g l = g' l (T,uj), where r = J \(t)dt and uj = z 2 , such 
that f = gl and g\ - g^ = 0. Therefore, 

u = \^ l {gl ] (T ) uj) + m\ ■ x)h % — \~ l (kt • x)k, 

p = 2A~ 2 (r7 • k t )g l (r,uj) + ±A" 2 (4 ■ k - 2k t ■ h)u 2 - (5.1) 

-\X~ l (n % ■ x)(fh\ t ■ x) — \\~ 2 (k ■ fh\ t )(h % ■ x)(k ■ x), 

where fh} • rh 2 — rh 1 ■ rh 2 = 0, k = fh 1 x rh 2 , n 1 = rh 2 X k, ft 2 = k x fh 1 , A = \k\ 2 , uj = k ■ x, 
t = J \{t)dt, and g\ - gi u = 0. 

For example, if fh = {r] l {t), 0, 0) and ft = (0,r/ 2 (t),0) with r] l (t) 7^ 0, it follows that 
u 1 = i^rHf 1 + r,] Xl ), u 2 = ( V 2 )-\f + r, 2 x 2 ), 



+ 



(UvWhiv'vT 1 - (faVM^Yr 1 ) 2 )*!, 



where f l = P(t,uj), f\ — f^ u = 0, r = J ' (r] l r] 2 ) 2 dt, and u = rfrfx?,. If fh 1 = (r/ 1 (t), J] 2 (t), 0) and 
rh 2 = (0,0,ry 3 (t)) with rj 3 (t) / and rf(t)rf{t) / 0, we obtain that 



u 1 



(77*77* ) — 1 77 1 (ffo; + rft x i) ~ V 2 ( r lt( r l 3 )~ 2uJ + tft x \ - ^i 3 ^)}, 
u 2 = (rfrf) -1 ^ 2 ^ + rf t Xi) + r] l (r]f(r] 3 )^ 2 uj + r} 2 x\ - r)\x 2 )^, 
u 3 = (r/ 3 )" 1 (/ + r/fx 3 ), 



p = 2(r] 3 ) 1 (r] 1 r] 2 — rjjr] 2 )(ri i r] t ) 2 g+^X 1 - 



•{A HivW - ZritVtWrf - 2r] 3 T] 3 rfrf t - 2(r l 3 ) 2 r? t rf t )J 2 + 
+ (r/ 3 ) 2 ((r/ 2 r/ 2 - r) X r)\ t )(x\ - x\) - 2{r, l tt rj 2 + r, l r, 2 tt )xix 2 )- 



vtt4}- 



Here / = f(r,u), f T - f WUJ = 0, g = g{r,uj), g T - g^ = 0, r = f(r) s ) 2 rfr) l dt, uj = ^(rfxi - r/ 1 ^), 
and A = (rf) 2 rfrf. 
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Note 5.1 The equation 

■ rh 2 — rh 1 ■ fa 2 = (5-2) 



can easily be solved in the following way: Let us fix arbitrary smooth vector-functions 
m 1 ,fe C°°((£o,£i),R 3 ) such that rh l (t) ^ 0, l(t) / 0, and rh l (t) • T(t) = for all t G (t Q ,tx). 
Then the vector- function fh 2 = rh 2 (t) is taken in the form 

rh 2 (t) = p(t)rh l + l(t). (5.3) 

Equation f |5.|jj implies 

p(t) = f(m l ■ m 1 )" 1 ^ 1 • f- m 1 • £)d£. (5.4) 
B. Let 0. By means of the transformation m 1 — ► a^fn? , where a,j = const and detjajj} = 



C, we make C = 1. Then we obtain the following solution of the NSEs (1.1) 
u= X' l (e ij {t)gl(T,uj) +6 i0 (t)io + mi-x - X^^kXm^-xfjn 1 - 

~\^ 1 (k t ■ x)k, 

p = 2\- 2 {n i ■ fctX^yfofcO + |^°(i)^ 2 )+ ( 5 - 5 ) 

+ ^A~ 2 (4 -k-2k t - k t )u? - \X- X {^ ■ x){m\ t ■ x)- 

— \\~ 2 {k ■ m\ t )(n l ■ x)(k • x). 

Here m\ ■ fh 2 — m 1 ■ m 2 = 1, k = rh 1 x rh 2 , n 1 = fh 2 x k, n 2 = k x m , A = \k\ 2 , uj = k ■ x, 
t = f X(t)dt, and g\ — g l UUJ = 0. (8 ll (t), 6 (t)) (i = 1, 2) are linearly independent solutions of the 
system 

Q\ + \- l (fn l ■ m 2 )6 l - A _1 (m* • rh l )9 2 = 0, (5.6) 

and (9 10 (t) , 6 20 (t)) is a particular solution of the nonhomogeneous system 

e\ + A _1 (m* • m 2 )6 l - A _1 (m* • m 1 )^ 2 = 2A~ 2 ((A? x • m l ). (5.7) 

For example, if fh 1 = (rj cos ip,rj sin tp,Q) and m 2 = (—n sin ip , rj cos , 0) , where 77 = r/(i) 7^ 
and ip = — I / (n)~ 2 dt (therefore, • fh 2 — rh 1 ■ fh\ = 1), we obtain 

u 1 = r? _1 (/ 1 cos ip — f 2 simp + n t x\ — \i]~ l X2), 



u: 



' ~ V simp + f 2 cos ip + r] t X2 + hv 



u 3 = -2r] t rj 1 x 3 , 
Here f = f(r,u;), f T - fL = 0, r = f(r/) 4 dt, and u> = ( V ) 2 x 3 . 
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Note 5.2 As in the case C = 0, the solutions of the equation 
ml ■ fh 2 — fh 1 ■ fh 2 = 1 



can be sought in form ( |5.^ ), As a result we obtain that 
p(t) = J Im 1 !" 2 ^ 1 • f- rh 1 ■ l t - l)dt. 



(5.8) 
(5.9) 



Note 5.3 System ^5.b\ ) can be reduced to a second-order homogeneous differential equation either 



in v , i.e., 



(Alm 1 ! -2 ^ 1 ) + ({{m 1 ■ fh 2 )\fh 1 \- 2 ) t + l™ 1 ^)^ = 
(then 2 = \fh 1 \~ 2 {\6 1 + {fh 1 • fh 2 )6 1 ) ), or in 2 , i.e., 

(\\m 2 \- 2 ef) t + (-{{rh 1 ■ m 2 )\rh 2 \-\ + |m 2 |- 2 )0 2 = 



(5.10) 



(5.11) 



(then 6 1 = \rh 2 \ 2 {—X9 2 + {rh 1 ■ fh )6 )). Under the notation of Note 5A equation ( 5.1C ) has the 
form: 



{{I- l)6 1 ) t + (m 1 !" 2 ^ 1 • I- m 1 ■ fr al 



0. 



(5.12) 



The vector-functions fh 1 and I are chosen in such a way that one can find a fundamental set 
of solutions for equation ( \5.1^) . For example, let fh x fht ^ Vi € {to,t\). Let us introduce the 
notation fh := fh 1 and put I = r]{t)fh x fh t , where r\ G C°°{{to, t±), R), rj{t) ^ Vi € {tQ,t\). 
Then 

fh ■ I = 0, fht • I — fh ■ If = 0, fh 2 = —{f \fh\~ 2 dt)fh + rjfh x fht, 
k = r]fh x {fh x fht), A = (7?) 2 |m| 2 |m x fht\~ 2 , 
h 2 = r]\fh\ 2 fh x fht, h 1 = {J \fh\~ 2 dt)h 2 + {r]) 2 \fh x fht\~ 2 fh, 
9 u {t) = J{r])- 2 \fh x m t \- 2 dt, 9 21 {t) = 1 - 9 U f \fa\- 2 dt, 

e 12 {t) = i, e 22 {t) = -f\fh\- 2 dt, 

9 10 {t) = 2 f{{{fh x fh t ) ■ fhtt)\fh x fh t \~ 2 + f ri~ 1 \fh\~ 4: dt)- 
• rr 2 \fh x fht\~ 2 dt, 

e 20 {t) = -e 10 {t) f \m\- 2 dt + 2 / n'^fh^dt. 

Consider the following cases: fhxfht = 0, i.e., fh = x(i)a, where x{t) G C°°{{to, ti),R), x{t) 7^ 
Vt € {to,t\), a = const, and \a\ = 1. Let us put 

l{t)=r ] 1 {t)b + r ] 2 {t)c, 

where rj X ,if £ C°°{{to, ti),R), {r, 1 (t) , rf {t)) ^ (0,0) Mt £ (i ,ii), b = const, \b\ = 1, a • b = 0, and 
c = a x 6. ITien 



-(x/x 2 dt)a + rr x b + rfc, k = x^c - xv 2 b, 
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A = (x) 2r ?*??% n 2 = (x) 2 (Vfr + r ? 2 c) ) n 1 = (/ X 2 dt)n 2 + xV l V l d, 
0™ = 2j(r^ri 1 -ri i 4)x~ 1 {v i rf)~ 1 dt, 9 20 = -9 10 J X ~ 2 dt. 



Note 5.4 In formulas (5A) and (5.1) solutions of the NSEs ([l.\ ) are expressed in terms of solu- 
tions of the decomposed system of two linear one- dimensional heat equations (LOHEs) that have 
the form: 



9t = 9L 



(5.13) 



The Lie symmetry of the LOHE are known. Large sets of its exact solutions were constructed [|/. 



The Q-conditional symmetries of LOHE were investigated in [H]. Moreover, being decomposed 
system ( 5. IS ) admits transformations of the form 



~g l {r',J) =F 1 (r,u;,g 1 (T,u;)), r' = G 1 (t,lu), J = H\t,u>), 

~g 2 (r\J') = F 2 (r,u;,g 2 (T,iv)), r" = G 2 (r,u), u» = H 2 (t,u), 

where (G 1 ,!! 1 ) ^ {G 2 ,H 2 ), i.e. the independent variables can be transformed in the functions g l 
and g 2 in different ways. A similar statement is true for system ( 5.1f\ )-( 5.2(\ ) (see below) if e = 0. 



Note 5.5 It can be proved that an arbitrary Navier- Stokes field (u,p), where 

U = W(t,ij) + (jfe*(t) • X)l l (t) 

with k\l % G C°°((t ,ii),K 3 ), k 1 x k 2 ^ 0, and u) = (k 1 x k 2 ) ■ x, is equivalent to either a solution 



from family ( \5.1\ ) or a solution from family (5.5). The equivalence transformation is generated by 
R(m) and Z(x)- 



5.2 Investigation of system (|37l3|) — (|37l^) 

Consider system 8 from Subsec. |3.2| , i.e., equations ( |3.13| )-( |3.16| ). Equation ( 3.1 6| ) immediately 



gives 



w 



-\ptP 1 + {v-^)z 2 2 , 



(5.14) 



where rj = rj(t) is an arbitrary smooth function of z\ = t. Substituting ( 5.14| ) into remaining 
equations (|5.13j )- (|5.15|) , we get 



+ (W 2 - X ) 2 Z2 3 , 
w l - w 22 + i 7 !^ 1 ~ \PtP~ X Z 2 )w\ = 0, 

w l ~ w 22 + (vz-i 1 ~ \ptP~ X z 2 )w\ + e(w 2 - x)z2 2 = 0. 



(5.15) 

(5.16) 
(5.17) 
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Recall that p = p(t) and x = x(t) are arbitrary smooth functions of t; e S {0; 1}. After the 
change of the independent variables 



r = J\p(t)\dt, z = \p(t)\^z 2 
in equations ( |5.16| ) and ( |5,17D , we obtain a linear system of a simpler form: 

2 w 2 zz + fj(T)z~ l w 2 z = 0, 



0, 



where ?)(r) = rj(t) and x(r) = x(t)- Equation ( |5.15| ) implies 



\{{ptp l )t~\{ptP l ) 2 )zl - r] t hx\z 2 \ 



(r, - l) 2 z 2 " 2 + JV(r, z) - x(r)fz^dz 2 . 



(5.18) 

(5.19) 
(5.20) 

(5.21) 



Formulas ( 5.14[) , (|5.18|) - (|5.21|) , and ansatz (3J3) determine a solution of the NSEs (1.1). 

If e = system ( |5.19| )-( [5.20| ) is decomposed and consists of two translational linear equations 
of the general form 



f T + fj(r)z- 1 f z -f z 



0.. 



(5.22) 



where fj = fj (fj = fj - 2) for equation (fyj|) ( (HQ )). Tilde over r\ is omitted below. Let us 
investigate symmetry properties of equation ( |5.22|) and construct some of its exact solutions. 

Theorem 5.1 The MIA of H5.2%) is given by the following algebras 

a) ^ =< fd f , g(T,z)d f > if t?(t) ^ const; 

b) L 2 =< d T , D, n, fd f , g(r, z)d f > if t](t) = const, rj {0; -2}; 

c) l 3 =< d T , d, n, d z + \i}z- l fd f , g = 2rd T -(z- nz-^fdf, fd f , 

g{r,z)d f > if r ? G{0;-2}. 

Here D = 2rd T + zd z , II = 4t 2 <9 t + Arzd z — (z 2 + 2(1 — ry)r)/9j; g = g(r, z) is an arbitrary solution 
of 

When rj = 0, equation ( |5.22j ) is the heat equation, and, when r\ = —2, it is reduced to the heat 
equation by means of the change f = zf. 

For the case rj = const equation (5.22) can be reduced by inequivalent one-dimensional subal- 
gebras of L 2 . We construct the following solutions: 

For the subalgebra < d T + afdj >, where a 6 {—1; 0; 1}, it follows that 

f = e - T z u {C l J u {z) + C 2 Y v {z)) if o = -l, 

f = e T z u {C l I u {z) + C 2 K u {z)) if a = l, 

f = C 1 z ri+1 +C 2 if a = and rj ^ -1, 

f = Ci\az + C 2 if a = and rj = — 1. 

Here J u and Y v are the Bessel functions of a real variable, whereas l v and if„ are the Bessel 
functions of an imaginary variable, and v = + 1). 
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For the subalgebra < D + 2afdf >, where a G E, it follows that 
/ = |r|° e -2>| ^-1)^(1^ _ l) _ a , l(r) + l),u) 
with oj = |z 2 t _1 . Here W(x, /j,u>) is the general solution of the Whittaker equation 
4w 2 W^ = (uj 2 - Akuj + 4/x 2 - 1)W. 

For the subalgebra < d T + U + afdf >, where a € K, it follows that 
/ = (4r 2 + ljsfa- 1 ) exp(-ro; + |a arctan 2r)v?(u;) 
with u = z 2 (At 2 + The function 92 is a solution of the equation 

Awip W u + 2(1 - 77)^ + - a)(/j = 0. 



For example if a = 0, then y?(o;) = u/^ (Ci J m (t;lj) + C2Y" M (^)j, where ji = |(ry + 1) 
Consider equation ( 5.22| ), where rj is an arbitrary smooth function of r. 



Theorem 5.2 Equation $5.22\ ) is Q- conditional invariant under the operat 

Q 1 =8 T + g\r, z)8 z + (g 2 (r, z)f + g 3 (r, z))d f 
if and only if 

g\ - r)z~ l g\ + r}z~ 2 g x - g\ z + 2g 1 z g 1 - r] T z~ l + 2g 2 = 0, 
g k T + nz^g k - g k zz + 2g\g k = 0, k = 2, 3, 

and 

Q 2 = d z + B(T,z,f)d f 
if and only if 



ors 



B T - r)z- 2 B + r,z- l B z - B z 



2BB 



zf 



B B 



if 



0. 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



An arbitrary operator of Q-conditional symmetry of equation (5.2i) is equivalent to either an op- 



erator of form (5.25) or an operator of form (5.25). 



Theorem 5.2 is proved by means of the method described in p3[. 



Note 5.6 It can be shown (in a way analogous to one in \lS\j) that system ( j5. 2 A) is reduced to the 
decomposed linear system 



r T+r] z-'r z -fi 







(5.27) 



by means of the following non-local transformation 

fzzf / fz 



9 



flf 2 - /7I 
f 1 f 2 - f 1 !' 2 

J zzJ z J z J zz 



:z +r ] z-\ 



(5.28) 



flf 2 - /7f ' 

/3 z -^- 1 /, 3 +5 1 /l-5 2 / 3 . 
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Equation ( 5.26 ) is reduced, by means of the change 

B = -* T /*f, * = *(r,z,/) 
and the hodograph transformation 

Vo =t, yi = z, y2 = $, * = /, 
to the following equation in the function Vl/ = ^(2/0,2/1,2/2)' 

Therefore, unlike Lie symmetries Q-conditional symmetries of ( |5.22| ) are more extended for 
an arbitrary smooth function n = T]{t). Thus, Theorem |5.2| implies that equation ( 5.22j ) is Q- 
conditional invariant under the operators 

8 Z , X = d T + (r ] -l)z- 1 d z , G=(2T + C)d z -zfd f 

with C = const. Reducing equation ( ]5.22 ) by means of the operator G, we obtain the following 
solution: 



f = C 2 (z 2 -2f(n(r)-l)dr)+C 1 . 

In generalizing this we can construct solutions of the form 

N 

f = Y / T k (r)z 2k , 

where the coefficients T k = T k (r) (k = 0,N) satisfy the system of ODEs: 
T T fc + (2k + 2)(t](t) -2k- \)T k+1 = 0, 

k = 0,N- 1, T? = 0. 



(5.29) 



(5.30) 



(5.31) 



Equation (5.31) is easily integrated for arbitrary N G N. For example if JV = 2, it follows that 
/ = C 3 {z 4 - 4* 2 /fa(r) - 3)dr+ 8/((7/(r) - 1) Mr) ~ 3)dr)dr} + 

+C 2 {z 2 -2f(n(r)-l)dT}+C 1 . 

An explicit form for solution (|5.30|) with N = 1 is given by ( 5.29f) . 
Generalizing the solution 



/ = C exp{-z 2 (4r + 2C)' 1 + J( v (t) - l)(2r + C^dr} 



(5.32) 



obtained by means of reduction of ( |5.22| ) by the operator G, we can construct solutions of the 
general form 



N 

f= ^(tXz&t + C)- 1 ) 2 *- 

■ exp{-z 2 (4r + 2C)- 1 + /^(r) - l)(2r + C^dr}, 
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(5.33) 



where the coefficients S k = S k (r) (k = 0,N) satisfy the system of ODEs: 
S k + (2k + 2)(t?(t) -2k- l)(2r + C)- 2 S k+1 = 0, 

k = 0,N- 1, = 0. 
For example if N = 1, then 

/ = {d(z 2 (2r + C)- 2 - 2J( V (t) - l)(2r + C)- 2 dr) + C }- 

• exp{-z 2 (4r + 2C)' 1 + /(ij(t) - l)(2r + C)- 1 ^}. 

Here we do not present results for arbitrary TV as they are very cumbersome. 

Putting g 2 = g 3 = in system ( 5.24j ), we obtain one equation in the function g 1 : 



9 l T ~ V z 1 9 1 z + V? V ~ glz + 2 ^5 fl ~VtZ 1 = 0. 
It follows that g 1 = —g z /g + (t) — 1)/^, where g = g(r, z) is a solution of the equation 

9t + (v~ ?)z~ l gz ~ 9zz = 0. 
Q-conditional symmetry of ( |5.22| ) under the operator 

Q = d T + (-g z /g+(v-l)/z)d z 
gives rise to the following 



Theorem 5.3 If g is a solution of equation ( 5.35 ) and 



f(r,z)= J z z Q z>g(T,z>)dz>+ 

+ i: {zo9z{r', zo) ~ (r?(r') - l)g(r>, z ))dr', 
where (tq,Zq) is a fixed point, then f is a solution of equation ( 5.2% ). 

Proof. Equation ( p. 35 ) implies 

(zg) T = (zg z - (77 - l)g) z 
Therefore, f z = zg, f T = zg z - (77 - l)g and 

f T + rjz~ 1 f z - f zz = zg z - (77 - l)g + ng - (zg) z = 0. QED. 
The converse of Theorem |5.3] is the following obvious 



Theorem 5.4 If f is a solution of ( |5. 22] ), the function 

9 = z~ X fz 
satisfies ( 5.3$) . 
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Theorems and 5A imply that, when r\ = 2n (n E Z), solutions of (|5.22|) can be constructed 



from known solutions of the heat equation by means of applying either formula ( |5.37| ) (for n > 0) 
or formula fl5,38| ) (for n < 0) |n| times. 

Let us investigate symmetry properties and construct some exact solutions of system ( |5,19| ) 
( 5.20| ) for e = 1, i.e., the system 



w\ - w\ z + f]{r)z- x w\ = 0, (5.39) 
w 2 - w 2 zz + (f)(r) - 2)z~ l w 2 z + (w 1 - x(r))z- 2 = 0. (5.40) 



If (w 1 ^ 2 ) is a solution of system ( |5.39| ) — ( 5.4C| ) , then {w l ,w 2 + g) (where g = g(r,z)) is also a 



solution of ( |5.39| )-( |5.40| ) if and only if the function g satisfies the following equation 

9r - 9zz + (r?(r) - 2)z~ 1 g z = (5.41) 
System ( 5.39| )-( [5.4C| ), for some x = x( t )j nas particular solutions of the form 



N N-l 

w ± = T k {T)z 2 \ w 2 =Y, S k {T)z 2 \ 

k=0 k=0 

where T°(t) = x(t). For example, if x(t) = ~~ 2C\ J(t}(t) — l)dr + C 2 and N = 1, then 
w 1 = d(z 2 - 2/(r?(r) - l)dr) + C 2 , w 2 = -dr. 
Let x(t) = 0. 



Theorem 5.5 The MIA of system fi5^ )- $54(\ ) with £(r) = is given by the following algebras 



a) < w l d w i, u> J (r, z)d w % > if fj(r) 7^ const; 

b) < 2rd T + zd z , d T , w l d w i, w 1 {t, z)d w i > if t/(t) = const, fj 7^ 0; 

c) < 2rd T + zd z , d T , w 1 z~ 1 d w 2, w i d w i, w 1 (t, z)d w % > if fj = 0. 
Here (w 1 ,w 2 ) is an arbitrary solution of ( |5. 39[ )- fl5.40l ) with x( T ) = 0- 



For the case x( T ) = and t}(t) = const system ( 5 . 39| ) — ( |5 . 40| ) can be reduced by inequivalent 



one-dimensional subalgebras of its MIA. We obtain the following solutions: 
For the subalgebra < d T > it follows that 
w 1 = C\ lnz + C 2 , 

w 2 = iCi(ln 2 z - hxz) + \C 2 In z + C 3 z~ 2 + C 4 

if fj = -i; 

w 1 = C 1 z 2 + C 2 , 

w 2 = \Ciz 2 + \C 2 In 2 z + C 3 In z + C 4 
if 77 = 1; 



w 



dz^ 1 + C 2 , 



w 2 = \C x {f) + l)" 1 ^ 1 + C 2 {f] - l)- 1 In 2 + C^~ x + C A 



if r? £{-!;!}. 
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For the subalgebra < d T — w l d w i > it follows that 
w 1 = e^z^^^iz), w 2 = e- T z^-i)^( z ) f 
where the functions ip 1 and ip 2 satisfy the system 

z 2 V4 + ziPl + {z 2 -\{f 1 + l) 2 )^ 1 = 0, (5.42) 
z 2 ip 2 zz + zip 2 + {z 2 - ±(r) - l) 2 )ip 2 = zip 1 . (5.43) 

The general solution of system ( |5.42| )-( |5,43| ) can be expressed by quadratures in terms of the Bessel 
functions of a real variable J u (z) and Y u (z): 

tP 1 = C 1 J v+1 {z) + C 2 Y 1/+1 (z), 

tP 2 = C 3 J v (z) + C 4 Y v (z)+ 

+%Y u (z) J J u {z)^ l {z)dz - \J v (z) J Y v {z)iP l {z)dz 

with v = | (77 — 1); 

For the subalgebra < d T + w l d w i > it follows that 
w i = e T z^ f>+1 >i{f L (z), w 2 = e T z^ f >- 1 ^ 2 (z), 
where the functions ip 1 and ip 2 satisfy the system 

z 2 *Pl z + z^ z -(z 2 + \{f l + l) 2 )^ 1 = 0, (5.44) 

z 2 ip 2 zz + zip 2 -(z 2 + l(fi- l) 2 )ip 2 = zip 1 . (5.45) 

The general solution of system ( |5.44j )-( |5.45| ) can be expressed by quadratures in terms of the Bessel 
functions of an imaginary variable I u {z) and K v (z): 

1P 1 = C l I u+l {z) + C 2 K u+l (z), 

1P 2 = C 3 I u (z) + C 4 K u (z) + 

K u (z) J I u (z)ip 1 (z)dz - I u (z)jK u (z)iP 1 (z)dz 

with v = 5(77 — 1). 

For the subalgebra < 2rd T + zd z + aw % d wl > it follows that 
w 1 = ItIV-^M^-^M. w 2 = \r\ a e-^\Lo\^- 3) tp 2 (Lo) 
with uj = |z 2 r _1 , where the functions ip 1 and ip 2 satisfy the system 

4u, 2 VL = + (« " J(V ~ 1)V + 1(V + I) 2 " lV, (5.46) 



4 W 2 VL = + (o - 1(V ~ 3))w + i(r} - l) 2 - l) V 
+2|w| 1 /2^i. 



2 + 



(5.47) 



The general solution of system ( 5.46| )- (|5.47| ) can be expressed by quadratures in terms of the 
Whittaker functions. 
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6 Symmetry properties and exact solutions of system 



As was mentioned in Sec. ||, ansatzes (|3,4| )-( ^7| ) reduce the NSEs ([O]) to the systems of PDEs of 
a similar structure that have the general form (see ( |3.12D ): 



WW; 



wh + si + a 2 w 2 



W l wf — W% + S2 — 012W 1 + a\w 3 



w i wf-wf i + a^w 3 + a 5 



(6.1) 



w 



«3, 



where a n (n = 1, 5) are real parameters. 

Setting ak = (k = 2, 5) in Q6.1|), we obtain equations describing a plane convective flow that is 



brought about by nonhomogeneous heating of boudaries [25]. In this case w l are the coordinates 
of the flow velocity vector, w 3 is the flow temperature, s is the pressure, the Grasshoff number A is 
equal to —a±, and the Prandtl number a is equal to 1. Some similarity solutions of these equations 
were constructed in [22]. The particular case of system ( |6.lD for a\ = ct2 = «4 = «s = and 
a3 = 1 was considered in plfl. 



In this section we study symmetry properties of system ( |6.1| ) and construct large sets of its exact 
solutions. 



Theorem 6.1 The 

1, 

2 
3 

4 

5 
6. 
7. 



Here D = Zidi 



Ex 


=< 3i, 


2 


E<2 






E3 


=<dt, 


o 2 


Ei 


=<dt, 


d 2 


E 5 


=<d u 


d 2 


Eq 


=<3l, 


2 


E~ 


=<d u 


d 2 


E% 


=< dx, 


d 2 



MIA of (6A) is the algebra 
d s > if a\ 7^ 0, 04 7^ 0. 

d s , d w s - aiz 2 d s > if a\ 7^ 0, a 4 = 0, (ai, a 2 , as) ^ (0, 0, 0). 
d s , d w 3 — a\Z2d s , D — 3w 3 d w 3 > if a\ / 0, a\~ = 0, k = 2, 5. 

d s , J, (w 3 + 05/04)5^,3 > if ot\ = 0, Q4 7^ 0. 

d s , J, d w 3 > if ai = a 4 = 0, (a 2 ,a 3 ) / (0,0), a 5 / 0. 

d s , J, d w 3, w 3 d w 3 > if ai = a 4 = a 5 = 0, (a 2 ,a 3 ) 7^ (0,0). 

d s , J, d w 3, D + 2w 3 d w 3 > if 7^ 0, ai = 0, I = 1,4. 

d s , J, d w 3, D,w 3 d w 3 > if a n = 0, n = 1,5. 

- w l d w i - 2sd s , J = zi&2 - z 2 d\ + w l d w 2 - w 2 d w i, <9; = d Zi . 



Note 6.1 The bases of the algebras Eq and Eq contain the operator w 3 d w 3 that is not induced by 
elements of A(NS). 

Note 6.2 If 7^ 0, the constant 0.5 can be made to vanish by means of local transformation 

w 3 = w 3 + 0.5/0.4, s = s — aia^a^ 1 Z2, (6-2) 

where the independent variables and the functions w % are not transformed. Therefore, we consider 
below that a§ = if 04 7^ . 
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Note 6.3 Making the non-local transformation 

s = s + a 2 ^, (6.3) 
where = w 2 , ^ 2 = —w 1 (such a function \& exists in view of the last equation of jjfi. 1\) ), in 



system (6.1 ) with a 3 = 0, we obtain a system of form ([6. 1\) with a 3 = a 2 = 0. In some cases 
(a\ / 0, 03 = 04 = «5 = 0, a 2 7^ 0; ot\ = a 3 = 04 = 0, a 2 7^ 0) transformation t \6.3( ) allows the 
symmetry of (\6. i|j to be extended and non-Lie solutions to be constructed. Moreover, it means that 



in the cases listed above system (6.1) is invariant under the non-local transformation 
Zi = e e Zi, w % = e~ £ w\ w 3 = e Se w 3 , s = e~ 2e s + a 2 (e~ 2e - 1)^, 

where 5 = —3 if 03 = a.4 = 05 = 0, a±, a 2 / 0; 

5 = 2 if a\ = a 3 = Q4 = 0, a 2 , 05 7^ 0; 
5 = if a\ = a 3 = 04 = 05 = 0, a 2 7^ 0. 

Let us consider an ansatz of the form: 
w 1 = a\(p l - a 2 tp 3 + biuj 2 , 

w 2 = a 2 (p l + anp 3 + b 2 0J 2 , 

w 3 = ip 2 + 63 U>2, 

s = h + d\0J 2 + d 2 ujiuj 2 + \d 3 uj 2 , 

where a 2 + a\ = 1, u = u\ = a\z 2 — a 2 z\, u> 2 = a\z\ + a 2 z 2 , B,b a ,d a = const, 

bi = Bai, b 3 (B + a 4 ) = 0, 

(6.5) 

d 2 = a 2 B - ai6 3 ai, d 3 = -B z - a\b 3 a 2) 



(6.4) 



Here and below (p a = f a (co) and h = h{uj). Indeed, formulas ( |6.4| ) and ( |6.5| ) determine a whole 
set of ansatzes for system ( |6.1[ ). This set contains both Lie ansatzes, constructed by means of 
subalgebras of the form 

< aid\ + a 2 d 2 + a 3 (d w a - aiz 2 d s ) + a 4 <9 s >, (6.6) 



and non-Lie ansatzes. Equation ( |6.5j ) is the necessary and sufficient condition to reduce ( |6.1| ) by 
means of an ansatz of form ( |6.3[) . As a result of reduction we obtain the following system of ODEs: 



f 3l fi ~ V>L, + MljV + di + d 2 w + a 2 ^ 3 = 0, 



</> 3< /4 - + hu>- a 2 Lp l + aian^ 2 = 0, 



V? 3 =<r, 



(6.7) 



where /in = —5, = —a\a 2 , fi 2 i = —b 3 , fi 22 = — «4, a = a 3 — B. If a = 0, system (6.7) implies 
that 

ip 3 = Cq = const, 
h = a 2 f (p {uj)duj — a\a\ f (p 2 (u>)dtu, 
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Table 1: Complete sets of inequivalent one-dimensional subalgebras of the algebras E\ — E% (a 
and ai (I = 1,4) are real constants) 



Algebra 


Subalgebras 


Values of 
parameters 


Ei 


< aidi + a 2 d 2 + a 3 d s >, < d s > 


a\ + a 2 = 1 


E2 


< aidi + a 2 d 2 + a 3 (d w 3 - aiz 2 d s ) >, 

< di + a 4 d s >, < d w s - a\z 2 d s >,< d s > 


a? + a 2 = 1, 
a 4 / 


E3 


< a\d\ + a 2 d 2 + a 3 (d w 3 - aiz 2 d s ) >, < d 1 + a 4 d s >, 

< D — 3w 3 d w 3 >, < d w z — a\z 2 d s >,< d s > 


a\ + a 2 = 1, 
a 3 G {-1;0; 1}, 
a 4 G {-1;1} 


E A 


< J + aid s + a 2 w 3 o w 3 >, < d 2 + a\d s + a 2 w°o w 3 >, 

< w 3 d w 3 + aid s >, < d s > 




E 5 


< J + aid s + a 2 o w 3 >, < d 2 + aid s + a 2 d w 3 >, 

< d w 3 + aid s >, < d s > 




Eq 


< J + aids + a 2 w 3 d w 3 >, < d 2 + aids + a 2 w 3 d w 3 >, 

< J + aids + a 3 d w 3 >, < d 2 + aid s + a 3 d w 3 >, 

< urd w 3 + aids >, < d w 3 + aids >, < d s > 


a 2 

a 3 G {-1;0;1} 


E 7 


< l) + aJ + 2w 3 d w 3 >,< J + aid s + a 2 d w 3 >, 

< d 2 + aids + a 2 d w 3 >, < d w 3 + a 2 d s >, < d s > 


a 2 G {-1;0; 1}, 
ai G {-1;0;1} 
ifa 2 = 


E$ 


< D + a J + a3W 3 d w 3 >, < -D + a J + a3d w 3 >, 

< J + aids + a4W 3 d w s >, < d 2 + aids + a 4 tt> 3 d^3 >, 

< J + aids + a 2 d w 3 >, < d 2 + aid s + a 2 d w 3 >, 

< w 3 d w s + aids >, < d w 3 + aids >, < d s > 


ai G {-1;0;1}, 
a 4 / 
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and the functions tp 1 satisfy system ( 4.23 ), where v\\ = d\ + 0:2 Co, ^21 = d 2 , ^12 = "5, ^22 = 0. If 
a 7^ 0, then v? 3 = era; (translating a;, the integration constant can be made to vanish), 

h = — \o 2 lo 2 + a 2 J tp l (uj)du) — a\a\ J (p 2 (ui)du!, 

and the functions satisfy system ( [4.29 ), where v\\ = d\, v 2 \ = d 2 + a 2 a, ^12 = «5, ^22 = 0. 

Note 6.4 Step-by-step reduction of the NSEs t \l-l\ ) by means of ansatzes (3.j i )-(3.7) and t \6.4j ) is 
equivalent to a particular case of immediate reduction of the NSEs (1A_) to ODEs by means of 



ansatzes 5 and 6 from Subsec. 1^.1 



Now let us choose such algebras, among the algebras from Table |l[ that can be used to reduce 



system ( |6.1[ ) and do not belong to the set of algebras (^6). By means of the chosen algebras we 
construct ansatzes that are tabulated in the form of Table 0. 

Substituting the ansatzes from Table ^ into system (|6.1| ), we obtain the reduced systems of 
ODEs in the functions ip a and h: 





- <pL, 










vl = 


0. 








- <pL> 




-</>L 




a 3 . 



ip 2 ip 2 -2h + aiip 3 sin a; + 2ip 2 



0.. 



+ h w — 2ip l w + a\(p 3 cos uj 

- 3</?V - 9^ 3 = 0, 

+ a 2 p 2 + K = 0, 

- 02V3 1 + 01 = 0, 

+ (a 2 <p 2 + a 4 - a 2 ,)^ 3 



0, 



(6.8) 



(6.9) 



o 11 

3. ujip tp 



a 2 a; V 2 + w 1 h UJ = 0, 

+ a 2 uj- 2 tp 2 - w~ + a 5 = 0, 
2V? 1 + wi^i = a 3 - 



(6.10) 



4. a;y>Vi 



-4,„2,„2 



a 2 uj~ 2 tp 2 + w _1 /i w = 0, 



<^V 2 



<#L + w V 2 - «2^ V + ai = 0, 



<^V 3 



+ a 2 a; 2 v?V 3 - w V 3 + («4 - ai w 2 )^ 3 
2V? 1 + wi^i = a 3 - 



(6.11) 
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Table 2: Ansatzes reducing system ( |6J| ) (r = (zf + z 2 ) 1 ^ 2 ) 



N 


Values 
of a n 


Algebra 


Invariant 
variable 


Ansatz 


1 


ol\ / 0, 
a k = 0, 
k = 2,5 


< £> - 3w 3 d w 3 > 


lo = arctan— 

n 


u; 1 = r _2 (z 1 V3 1 - z 2 <^ 2 ), 
w 2 = r ~ 2 (z 2 ip 1 + <p 2 ), 


2 


a\ = 0, 
a 5 = 


< d 2 + aid s + a,2W 3 d w 3 >, 
a 2 #0 


UJ = Z\ 


11 2 2 

W = if , w = <p , 
w 3 = ip 3 e a ' 2Z ' 2 , 
s = h + aiZ2 


3 


ai = 0, 
«4 = 


< J + aid s + a 2 <9 M ,3 > 


ijj = r 


1 1 —2 2 

wr = zi</? — z 2 r (p , 
w 2 = Z2f 1 + zir~ 2 ip 2 , 
■w 3 = ip 3 + a 2 arctan 
s = /i + ai arctan— 


4 


ai = 0, 
a 5 = 


< J + aids + a2W 3 d w 3 > 
a 2 ^ if 04 = 


UJ = r 


1 1 —2 2 
10 = Zl<p — Z 2 r if , 

w 2 = + zir~ 2 if 2 , 

rt q 0,2 arctan— 

e — /1 1 /i-i q it*"! - n n 


5 


a 5 / 0, 

«i = o, 

I = 1,4 


< D + aJ + 2w 3 d w 3 > 


to = arctan— — 
—alnr 


= r _2 (zi(/p 1 - z 2 (p 2 ), 
w 2 = r _2 (z 2 ^ 1 + zi(p 2 ), 

'3 2 H —2 7 
uj — / \u . j — / it 


6 


a n = 0, 
n = 1, 5 


< D + a J + ai<9„,3 > 


to = arctan— — 
—alnr 


w 1 = ^(znp 1 - z 2 (f 2 ), 
w 2 = r~ 2 (z 2 <p l + zitp 2 ), 

= </7 3 + ai lnr, 
s = r~ 2 h 


7 


a„ = 0, 
n = 1, 5 


< L> + a J + aiw 3 d w 3 >, 
ai #0 


a; = arctan— — 
—a In r 


w 1 = r~ 2 (ziLp l - z 2 (p 2 ), 
w 2 = r ~ 2 (z 2 ip l + zitp 2 ), 
w 3 = r ai (p 3 , s = r~ 2 h 
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5. 


(v 2 


- atp^ipi - 


(1+a 2 )^ 


- <pV - 


- 93 2 v? 2 — ah w — 2h = 


0, 




& 2 


- a^ifl - 


(i + « 2 )^L 


- 2(a^ 2 


+ ipl) + K = 0, 






(<P 2 


- a^)vl ~ 


(i + « 2 VL 


+ 2pV 


- 4v9 3 + 4a</9 3 + a 5 = 


= 






- ayl = 0. 










6. 


^ 


- a^fl - 


(l+a 2 )^ 


- *V - 


- c^V 2 - aha, — 2h = 


0. 






- W 2 - 


(l+a 2 )^L 


" 2(a^ 2 


+ ^i) + K = o, 








- a^vl ~ 


(1 + a 2 )^ 


+ aiy? 1 = 


= o, 








- afl = 0. 










7. 


(<f 2 


- a<p l )<pl - 


(l+a 2 )c^ 


- *V - 


- ip 2 ip 2 - ah u -2h = 


o. 




(^ 2 


- a^)vi - 


(i + « 2 VL 


- 2(a^ 2 


+ ifi) + K = o, 






(^ 2 


- a^)^l ~ 


(i + « 2 )^L 


+ aiv?V 


3 - a 2 (/? 3 + 2aa 1 ifl = 


o, 




^ 2 - 


- aipl = 0. 











(6.12) 



(6.13) 



(6.14) 



Numeration of reduced systems fl6.8Q - (|6.14 ) corresponds to that of the ansatzes in Table ||. Let 



us integrate systems (|0|)-( |6"7l4| ) in such cases when it is possible. Below, in this section, = const 
(A; = 1^). 

1. We failed to integrate system ( |6.8| ) in the general case, but we managed to find the following 
particular solutions: 

a) 95 1 = -6p( W + C 3 ,i(4-2C 1 ),C 2 )-2, 
tf? = tp3 = 0, h = 2< / o 1 + Ci; 

b) cp 1 = -6C 2 e 2C ^p(e c ^ + C 3 , 0, C 2 ) + 3C 2 - 2, 
y, 2 = 5Ci, v? 3 = 0, 

/» = -12C 1 2 e 2C ^p(e c ^ + C 3 , 0, C 2 ) - 2 - f C 2 - |C X 4 ; 

c) ^ = Ci, </? 2 = C 2 , </? 3 = 0, fr = -f(C? + Cf). 

Here p(t,hx,K2) is the Weierstrass function that satisfies the equation (see fl9| ): 

{p T f =4p 3 -x 1 p-x 2 - (6.15) 



2. If a3 = 0, the last equation of ( |6.9| ) implies that (/j 1 =C\. It follows from the other equations 
of (0) that 

V? 2 = C 3 + C 2 e c ^ - (aiCf 1 - a 2 )a;, 
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h = C 6 - a 2 C 3 uj - a 2 C 2 C^ 1 e Clu ' + ±a 2 (aiCf 1 - a 2 )uj 2 



if Ci ^ 0, and 

ip 2 = C 3 + C 2 uj + iaiw 2 , 

h = Cq- a 2 C-iUJ - \a 2 C 2 uj 2 - ^a 2 aiuj 3 
if Ci = 0. The function tp s satisfies the equation 



2\,„3 



Cl<ftj + ( a 2 - "4 - «2<£> )<P 



0. 



(6.16) 



We solve equation (6.16) for the following cases: 
A. C 2 =a 1 - a 2 C± = 0: 



<P 



C'2 



(C^" u + C 5 e 



H > 0, 
/i = 0, 



E2 



(C 4 cos^-//) 1 ^) + C 5 sin((- / u) 1 / 2 w)), /z < 0, 



where /i = \C 2 — a 2 . + 04 + a 2 C3. 

B. d = 01 = 0, C 2 ± (@) : 

^ = e 1 /^i/3(|(-a 2 C 2 ) 1 /^3/2 ); 

where £; = w + (C3a 2 — a| — a^)j{a 2 C 2 ). Here Z v (t) is the general solution of the Bessel equation 

C. d = 0, ai ^0 (||): 

^ 3 = (w + C 2 ar 1 )" 1 / 2 Ty(i/, \, {\ ai a 2 y l l 2 {u + CW) 2 ), 

where v = 5(^aia2) _1 / 2 (a 2 — «4 — a2C*3 + ^^C^di 1 ). Here W(x,h,t) is the general solution of 
the Whittaker equation ( [4.21] ). 

D. d + 0, C 2 ^ 0, ai - a 2 d = 

p? = e\^Z v {2Cl\-a 2 C 2 fl 2 e\ c n, 
where v = (7-f 1 (C 1 2 + 4(a4 + a 2 C l 3 — a 2 )) 1 ^ 2 . Here Z v (t) is the general solution of the Bessel equation 

(ED- 

E. d / 0, 01 - a 2 Ci / 0, C 2 = (||): 

^3 = e \C^ e /2 Zi/ ^2 {a2{aiC -l _ a2 ))V2 e 3/2) ) 

where ^ = uj + (a 2 — \C 2 — C 3 a 2 — Q4) / '(a 2 (aiC-f 1 — a 2 )). Here Z v {r) is the general solution of the 
Bessel equation ( 4.22 ). 

If a 3 0, then ip 1 = a 3 uj (translating u>, the integration constant can be made to vanish), 
<p 2 = C 1 + C 2 J e^ 2 doj + 01 J e^ 2 e -\<*^ 'dujdu + a 2 u, 
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h = C 3 - + aj)oj 2 - a 2 Ciuj - a 2 C 2 (u J e^ aW ' dw - 1 e^ a3t 



a 2 ai\ 



and the function ip 3 satisfies the equation 

vlu ~ a 3Utpl + (a 2 , - a 4 - a 2 tp 2 )ip 3 = 0. 
We managed to find a solution of ( |6.17 ) only for the case a\ = C 2 = 0, i.e., 



(6.17) 



^3 = e i^V (a i/2 (w + 2a 2 a 2 af), v), 

where v = 4a 3 ~ 1 (a:4 + a 2 C\ — a 2 (a 2 a^ 2 + 1)). Here V(t,i/) is the general solution of the Weber 
equation 



4Kr = (r 2 + v)V. 
3. The general solution of system ( |6.10D has the form: 

^ = du- 2 + ia 3 , 

ip 2 = C 2 + C 3 /w Cl+1 e* asw2 dw - \a 2 u 2 + 
ip 3 = C 4 + C 5 Juj c ^- 1 e l i a3Uj2 dio+ 

h = C 6 - \aluj 2 - \C\uj~ 2 + /(v? 2 (a;)) 2 a; _3 cL; - a 2 / w _ V 2 M d w. 



(6.18) 
(6.19) 
(6.20) 



(6.21) 



4. System fl6.11| ) implies that the functions ip 1 and h are determined by (|6.19| )-( |6.21]) , and the 
function ip s satisfies the equation 

<fL,- ((C l -l)u- 1 + \a 3 u)vl + (a 2 u;- 2 (a 2 -ip 2 ) - a 4 )v? 3 = 0. (6.22) 

We managed to solve equation ( |6.22 ) in following cases: 
A. C 3 = ai = 0, a 3 / 0: 

^ = ui^e^Vfo/i, \a 3 u; 2 ), 

where x=\{2-C 1 - (4a 4 + 2a 2 a 2 )a 3 ~ 1 ), \i = ±(C 2 - 4a 2 , + 4a 2 C 2 ) 1/2 . Here PF(x,/x,r) is the 
general solution of the Whittaker equation ( |4.21 ) . 
Let a 3 = 0, then 

' C 2 + C 3 lncj + i(ai + 2a2)u; 2 , C x = -2, 

t^ 2 = J C 2 + iC 3 w 2 + 5aiw 2 (lnw-i), Ci = 0, 

, <?2 + C 3 (C! + 2)-!^+ 2 - \C^\ ai - a 2 C 1 )u 2 , d ^ 0, -2. 
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B. C 3 = a\ — a 2 C\ = 0: 



>^(C 5 u v + C 6 u;- v ), [jl = Q,v^0, (6.23) 

4 w5 c i(C 5 + C 6 lnw), ^ = 0,^ = 0, 

where n = — a 4 , i/ = ±(C 2 - 4a 2 , + 4a 2 C 2 ) 1/2 . Here and below Z„ (r) is the general solution of the 
Bessel equation (|4.22[) . 

C. C3 = 0, C\ 7^ 0: (/? 3 is determined by ( |6.23D , where 

\i = \a 2 C{ l (a 1 - a 2 C 1 ) - a 4 , v= \{C\ - 4a 2 , + 4a 2 C 2 ) 1 / 2 . 

D. C\ = a\ = 0: c^ 3 is determined by ( |6.23 ), where 
fJ> = —\a 2 Cz — a 4 , f = (— a| + a 2 C 2 ) 1/2 . 

E. C 3 / 0, Ci {0; -2}, a 2 (ai - a 3 Ci) - 2a 4 Ci = 0: 
<p 3 = ujh c ^Z u {fxuj 1+ ^ Cl ), 

where n = 2C 3 1/2 (d + 2)~ 3 / 2 , v = (Ci + 2)~ 1 (<?i - 4a 2 . + 402C2) 1 / 2 . 

F. Ci = -2, C 3 / 0, 02(01 + 2a 2 ) + 4a 4 = (|l9[): 

where £ = In a; + C^ 1 (a 2 — a 2 C 2 — 1). 

G. Ci = 2, C 3 < 0, 1 - a 2 . + a 2 C 2 > 0: 

^ 3 = ^,^i(-C 3 ) 1 / 2 W 2 ), 



where 



|(-C 3 )~ 1 / 2 (-4a 4 + a 2 - 2a 2 a 2 ), 



/i 



|(1 -a\ + a 2 C 2 ) l l 2 . Here Wfo/i.r) is the 



general solution of the Whittaker equation ( 4.21| ). 

5-7. Identical corollaries of system ( |6.12 ), ( 6.13; ), and ( 6.14| ) are the equations 



aip 1 + Ci, 



h = a(l + a 2 )^ + (2 + 2a 2 - aCi)^ 1 + C 2 , 

(1 + a 2 Vl + (4a-dVi + yV+ V+(l + a 2 )^(C 2 +2C 2 ) = 0. 



(6.24) 
(6.25) 
(6.26) 



We found the following solutions of (6.26): 
A. If (l + a 2 )- 1 ^! +2C 2 ) <4: 

<pi = (4 - (1 + a 2 )-\C 2 + 2C 2 )) 1/2 - 2. 



B. If Ci = 4o: 



, r 4 (C 2 + 2C 2 ) " 
(l + a 2)i/2 + ° 4 > 3 3(1 + a 2 ) ' 3 



(6.27) 



(6.28) 
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Here and below p(r, x%, ^2) is the Weierstrass function satisfying equation (|6.15| ). If C2 = 2 — 6a 2 
and C3 = 0, a particular case of ( 6.28 ) is the function 

ip 1 = -6(1 + a 2 )uj 2 - 2 (6.29) 

(the constant C4 is considered to vanish). 

C. If d + 4a, (1 + a 2 Y x (C\ + 2C 2 ) - 4 = -9fi 4 : 

p 1 = -6 / u 2 e" 25 p(e^ + C 4 , 0, C 3 ) + 3/i 2 - 2, (6.30) 

where $ = (1 + a 2 )" 1 / 2 ^, fj, = \{Aa - Ci)(l + a 2 )" 1 / 2 . If C 3 = 0, a paticular case of flOOl) is 
the function 

= _6^ 2 e - 2? (e~ 5 + C 4 )~ 2 + 3// 2 - 2, (6.31) 

where the constant C 4 is considered not to vanish. 

The function <p 3 has to be found for systems (6.12), ( 6.13j ), and ( |6.14 ) individually. 

5. The function ip 3 satisfy the equation 

(1 + a 2 )<p 3 ^ - (d + 4a)<pl - (2^ - 4)<^ 3 - a 5 = 0. 
If is determined by ( |6.27| ), we obtain 
<p 3 = exp(|(l + a 2 )" 1 ^ + 4a)w)- 

C 5 exp(z// 2 c,;) + C , 6 exp(-*// 2 w), zy > ] 
C 5 cos((-^) 1 /2 a; ) +C7 6 sin((-^) 1 /2 t(J ) ) „<o + 

C 5 + C 6 w, ^ = J 

-a^ip 1 -4)- 1 , 

-a 5 (4a + Ci)- 1 a;, 2^-4 = 0, Ci + 4a / 

^^(l + a 2 )- 1 ^ 2 , 2^-4 = 0, Ci + 4a = 

where z/ = \{1 + a 2 )" 2 (Ci + 4a) 2 - (1 + a 2 )" 1 ^ - 2^). 

6. In this case ip 3 satisfy the equation 
(1 + a 2 )^ - C x ipl = ai ip\ 

Therefore, 



f 3 = C 5 + C 6 exp((l + a 2 y 1 Ciuj) + a 1 C\ 1 (/ p 1 (u)duj+ 
exp((l + a 2 )-^!^) / exp(-(l + a 2 )- 1 ^)^ 1 ^)^) 



for Ci 7^ 0, and 

(/? 3 = C5 + C^w + ai(l + a 2 ) -1 ^ / </3 1 (a;)cia; - / u:p 1 (uj)dui) 
for Ci = 0. 

7. The function y? 3 satisfy the equation 

(1 + a 2 )^ - (Ci + 2aia)(^ 3 + (aj - a W 3 = 0. (6.32) 
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A. If p 1 is determined by ( |6.27|) , it follows that 
p 3 = exp(|(l + a 2 ) _1 (Ci + 2a 1 a)u)- 

C 5 exp(zW 2 u;) + C 6 exp(-z// 2 w), v > 
C 5 cos((-^) 1 /2 a; ) + C7 6 S in((-^) 1 / 2 w), v < 
C 5 + C 6 w, 1/ = 

where v = \{l + a 2 )~ 2 {Ci + 2a ia ) 2 - (1 + a 2 ) _1 (ai - ai^ 1 )- 

B. If C\ = 4a, that is, 99 1 is determined by fl6.27| ), we obtain 

93 3 = exp(a(ai + 2)(1 + a 2 ) _1 a;)6l(r), 

where r = (l + a 2 ) _1 / 2 a; + C4. Here the function 9 = 6{t) is the general solution of of the following 
Lame equation ( [p| ) : 

9 TT + (6aip(r) + a? + 2ai - a 2 (2 + ai) 2 (l + a 2 )" 1 )^ = 

with the Weierstrass function 

p(r) = p(r, 1(4 - (1 + a 2 )- 1 !^ 2 + 2C 2 )), C 3 ) . 

Consider the particular case when C 2 = 2 — 6a 2 and C3 = additionally, i.e., 99 1 can be given 
in form ( |6.29| ). Depending on the values of a and ax, we obtain the following expression for ip 3 : 

Case 1. ax / —2, ax 7^ 2a 2 : 

3 1/2 ( a(2 + a 1 ) \ ( {{2 + a 1 ){a 1 -2a 2 )) 1 ' 2 \ 

where = (| — 6a!) 1 / 2 . 

Case 2. ai = —2: y? 3 = C5U; 4 + Cqlo^ 3 . 

Case 3. ai = 2a 2 : 

Case 3.1. 48a 2 < 1 : p 3 = |w| 1 / 2 e 2aw (C 5 a; <T + C 6 uj~ a ), 
where a = |vl — 48a 2 . 

Case 3.2. 48a 2 = 1, that is, a = iy^v^: ip 3 = \uj\ 1/2 (C 5 + C 6 law). 

Case 3.3. 48a 2 > 1 : p 3 = |w| 1 / 2 e 2a ^(C 5 cos( 7 hxu) + C 6 sin( 7 lnw)), 
where 7 = l\/48a 2 — 1. 

C. Let the conditions 

Ci/4a, (l + a 2 )- 1 (C 1 2 + 2C 2 ) -4 = -9/ 



be satisfied, that is, let p 1 be determined by ( |6.30 ). Transforming the variables in equation (|6.32j ) 
by the formulas: 



^3 = T -x/2 exp (i (d + 2aax)(l + a 2 )~ l uj)e{T) 



t = exp(-/z(l + a 2 )~ 1/2 w), 
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we obtain the following equation in the function 8 = 8{t): 
t 2 6 tt + (6 ai T 2 p(T + C 4 , 0, C 3 ) + o)0 = 0, 



(6.33) 



where a = n~ 2 {a\ + 2a% - \{l + a 2 ) _1 (C 2 + 2aai) 2 ) - 3ai + |. If a = 0, equation ( |6.33|) is a Lame 
equation. 

In the particular case when ip 1 is determined by (6.31), equation (6.33) has the form: 



r 2 (r + C A ) 2 9 TT + (6air 2 + <t(t + C 4 ) 2 )# = 0. 
By means of the following transformation of variables: 



(6.34) 



where z^i(^i — 1) + a = and ^2(^2 — 1) + 6ai = 0, equation ( 6.34| ) is reduced to a hyper geometric 
equation of the form (see (n| ) : 

S(Z - + (2(^1 + ^ 2 )£ - 2i/i)^ + 2^i/ 2 ^ = 0. 
If cr = 0, equation ( |6.34[ ) implies that 

(r + C 4 ) 2 6» rT + 6ai0 = 0. 
Therefore, 

= C7 5 | T + C^ 2 " 1 ' + C 6 |r + C 4 | 1/2+l/ 
if a\ < , where v = (~ — 6a!) 1 / 2 , 

= |T + C7 4 | 1 / 2 (C 5 + C 6 ln|r + C 4 |) 
if ai = 24 , and 

= |r + C 4 | 1 / 2 (C7 5 cos(^ln|r + C 4 |) + C 6 sin(^ln|r + C 4 |)) 
if ai > ^ , where f = (6ai — |) 1//2 . 

7 Exact solutions of system ( |2.9|) 



Among the reduced systems from Sec. only particular cases of system (2^) have Lie symmetry 
operators that are not induced by elements from A(NS). Therefore, Lie reductions of the other 
systems from Sec. || give only solutions that can be obtained by means of reducing the NSEs with 
two- and three-dimensional subalgebras of A(NS). 

Here we consider system ( |2.9| ) with p 1 vanishing. As mentioned in Note |2.5| , in this case the 
vector-function rh has the form rh = rj(t)e, where e = const, |e| = 1, and rj = r/(t) = \rh(t)\ 7^ 0. 
Without loss of generality we can assume that e = (0,0, 1), i.e., 

m= (0,0,J7(t)). 

For such vector fh, conditions are satisfied by the following vector h % : 
fl 1 = (1,0,0), n 2 = (0,1,0). 
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Therefore, ansatz fl2.4|) and system ( |2.9| ) can be written, respectively, in the forms: 
u 1 =v 1 , u 2 = v 2 , u 3 = (r](t))~ (v 3 + r) t (t)x3), 

(7-1) 

p = q - htt(t)(v(t)) xl 

where u = u (2/1,2/2,2/3), 9 = 9(2/1,2/2,2/3), Vi = Xi, Vs = t , and 

*; 3 + - u^- = 0, (7.2) 
uj + P 3 = 0, 

where p 3 = p 3 (t) = r\tjr\. 

It was shown in Note |2,8| that there exists a local transformation which make p 3 vanish. There- 
fore, we can consider system (|7.2|) only with p 3 vanishing and extend the obtained results in the case 
p 3 7^ by means of transformation ( |2.12| ). However it will be sometimes convenient to investigate, 
at once, system ( [7.2p with an arbitrary function p 3 . 

The MIA of ( |T72|) with p 3 = is given by the algebra 

B =< R 3 (ij), Z\X), Dl d t , Jl 2 , d v3 , v 3 8 v3 > 

(see notations in Subsec. [2.1|) , We construct complete sets of inequivalent one-dimensional subal- 
gebras of B and choose such algebras, among these subalgebras, that can be used to reduce system 
( |7.2| ) and do not lie in the linear span of the operators R^ty), Z 1 (X), J\ 2 , i-e., the operators 
which are induced by operators from A(NS) for arbitrary p 3 . As a result we obtain the following 
algebras (more exactly, the following classes of algebras): 
The one-dimentional subalgebras: 

1. B\ =< D\ + 2xJ 1 1 2 + 2-/v 3 d v s + 2(3d v s >, where 7/3 = 0. 

2. B\ =< d t + *J\ 2 + -yv 3 d v 3 + j3d v a >, where 7/3 = 0, x G {0; 1}. 

3. B\ =< Jl 2 + jv 3 d v3 + Z 1 (X(t)) >, where 7 / 0, AG C°°((t , h), R). 

4. B\ =< R 3 $(t))+jv 3 d v 3 >, where 7 ^ 0, $(t) = ty\t),il?{t)) + (0,0) Vi G (i ,ii), 

^GC°°((i 0! ii),M). 
The two-dimentional subalgebras: 

1. B\ =< d t + (3 2 d v a, D\ + kJ\ 2 + 7^3 + fi x d v 3 >, where 7ft =0, (7 - 2)/3 2 = 0. 

2. fl| =< D\ + 2~f lV 3 d v 3 + 2/3ia„ 3 , J{ 2 + 7 2 w 3 <9„3 + /? 2 d„3 + >, where 

7i/3i = 0, 7 2 /3 2 = 0, 7i/3 2 - 72ft = 0. 

3. S| =< + 2xJjj + 27iu 3 ^3 +2/3i^s, Ifefltf+Va cost, \t\ (T+l l 2 sinr) + 72V 3 <9„3 + 

(3 2 d v 3 + Z 1 (e\t\ a ' 1 ) >, where r = xln|i|, (71 + a)/3i - 72 /?i = 0, 0-72 = 0, ea = 0. 

4. B| =< d t + 7iw 3 9 t ,3 + /3i<9 v s, j\ 2 + 7 2 w 3 3„3 + /3 2 <9„3 + Z 1 (e) >, where 71 ft = 0, 

= 0, 71/02 - 72#L = 0. 

5. S| =< fit + xJ| 2 + 7iu 3 a,,3 + A^s, R3(e at cosxt,e at sinxi) + Z^ee* 7 *) + r y 2 v 3 d v 3 + 

/3 2 <9„3 >, where (71 + a)(3\ - = 0, 0-72 = 0, ecr = 0. 
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6. B 2 =< R^ 1 ) +jv 3 d v s, R 3 (^ 2 ) >, where ^ = {^{t),^(t)) + (0,0) Vt G (t ,*l), 
^ij e ^((tcii),!!), ^ • V^ 2 - V^ 1 • = 0, 7^0. Hereafter ij; 1 ■ $ 2 := ^ u ^ 2i . 



Let us reduce system ( |7.2| ) to systems of PDEs in two independent variables. With the algebras 
B\-B\ we can construct the following complete set of Lie ansatzes of codimension 1 for system 
(p|) with p 3 = 0: 



1. v 1 = \t\ ^(u? 1 cost — u; 2 sinr) + iyii 1 — xy%t l , 
v 2 = |t| _1 / 2 (w 1 sinr + w 2 cost) + \y2t~ 1 + xyit^ 1 , 

(7.3) 

v 3 = \t\~<w 3 + (3 In \t\, 

q = |t|-l s+ l(^ + l) t -V, 
where t = >r In | £ | , 7/3 = 0, 

z\ = |t| _1 / 2 (yi cos t + y2 sinr), Z2 = |i| _1 ^ 2 (— yi sinr + 2/2 cos r). 
Here and below w a = w a (zi, Z2), s = s{z\, Z2), r = (y 2 + y 2 ) 1 ^ 2 . 

2. v 1 = w 1 cos xt — w 2 sin xt — xy2, 

v 2 = w 1 sin xt + w 2 cos xt + xy\ , 
v 3 _ ^3 e 7t + ^ 

q =s + i>f 2 r 2 , 
where ;nr G {0; 1}, 7/5 = 0, 

^1 = y\ cos >rf + 2/2 sin xt, Z2 = — yi sin xt + y 2 cos xtf. 

3. v l = y\r~ l w 3 — y2r~ 2 w l — r yy 2 r~ 2 , 



(7.4) 



v 2 = y2T l w 3 + y\r 2 w 1 +^y\r 2 , 

■y 3 = ^2 gTarctan 2/2/2/1 

g = s + A(t) arctan 1/2/2/1, 
where z x = t, z 2 = r, 77^0, A G C ,oo ((t ,ti) J M). 

4. « =(^-^)- 1 ((«; 1 +7)^ + «^e+(^-y)^ ( -«2^) 



(7.5) 



■u 3 = ui 2 exp(7(^ • V>) 1 (^-J/)) (7-6) 

<? = « - (-0 ■ ^)" 1 (V'tt • y)(^ • y) + 3$ • V0~ 2 (V^ ■ • y) 2 , 

where a* = t, z 2 = (9-y), 7 r 4 0, ^(w 1 ,*; 2 ), y = (yi,y 2 ), ft G C°°((io, K), = (-V 2 ,V' 1 )- 



Substituting ansatzes (|7.3| ) and (7^) into system ( |7T2; ) with p 3 = 0, we obtain a reduced system 



of the form (6.1), where 

a\ = 0, 012 = —1, 03 = — 2xr, 04 = 7, 05 = /3 if t > and 
«i = 0, 02 = 1, «3 = 2x, 04 = —7, «5 = — j5 if £ < 
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for ansatz ( |7.3|) and 

OL\ = 0, «2 = 0, 03 = — 2xr, 04 = 7, 05 = /3 



for ansatz (7/4). System ( |6.1|) is investigated in Sec. y in detail. 

Because the form of ansatzes ( |7.3| ) is not changed after transformation ( |2.12| ) , it is convinient to 
substitute their into a system of form ( |7.2[) with an arbitrary function p 3 . As a result of substituting, 
we obtain the following reduced systems: 



9 ^9 9 19 219 

W{ + W W2 — W22 — Z 2 W2+ 7-^2 W W 



W 2 + Z 2 w3 



3. wf + W 3 W2 — Z 2 3 (w 1 + 7) 2 — (^22 + z 2 lw 2 ~ z 2 2llj3 ) + s 2 = 0, 

w\ + w 3 Wr, — W22 + z 2 X w\ + A = 0, 

0, 

-771/77. 

4. u>{ + li; 3 ^ — (V> • ^)w22 = 0) 

U> 3 + U) 3 ^ — (tp ■ 1p)w22 + (V> ■ i>)s2 + 2(w 1 + 7)(V' • 0)(lp ■ 

2(& • ■ + (2& • & - ^ ■ ■ ^-^2 = 0, 

+ U) 3 7i>2 - • V>)u>22 + 7(^ • VO" 1 ^ 1 + ■ 0)(lp ■ tp)~ 1 Z2)w 2 

w 2 +Vt/v = 0- 



(7.7) 



(7- 



Unlike systems 8 and 9 from Subsec. ^2, systems (7/7) and ( JLaj ) are not reduced to linear systems 
of PDEs. 

Let us investigate system (|7.7j ). The last equation of ([7.7[) immediately gives 

(w| + z 2 " 1 ^ 3 )2 = w| 2 + ZoT 1 ^ - Z2~ 2 w 3 = 0, io 3 = (x - 1)^2 1 - iWT 1 -^, (79) 

where x = x(t) 1S an arbitrary differentiable function of t = Z2- Then it follows from the first 
equation of (|7.7|) that 



Z2\- 



J z 2 " V + 1 fdz 2 -\{ X - l?Z2 2 + \z\ ((mh)t - Hvt/V?) - Xt In 



Substituting ( [7.91 ) into the remaining equations of (7/7), we get 
w l ~ w h + {X z 2 l - \ritr]- 1 z 2 )w\ + A = 0, 



w l ~ w 22 + ((X ~ 2)^2 1 - htV L Z2)W2 + 7^2 



-2„.,1„.,2 



0. 



(7.10) 



By means of changing the independent variables 

T = J\rj(t)\dt, z = \i ] (t)\ 1/2 Z2, 
system (7.10) can be transformed to a system of a simpler form: 
wl ~ «4 + Xz~ l w 2 z + \\fj\~ 1 = 0, 



(7.11) 



w'j z + (x — 2)z 1 w'j + jz 2 w 1 w 2 



0: 



(7.12) 
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where 77 (r) = r](t), x(t) = x(*), and A(r) = X(t). 

If A(t) = — 2Cr](t)(x(t) — 1) f° r some fixed constant C, particular solutions of ( [7.10 ) are functions 

w 1 = Cr)(t)z%, w 2 = f(z!,z 2 ) exp(-fC Jrj(t)dt), 

where / is an arbitrary solution of the following equation 

h ~ /22 + ((x " 2)^ - htr j - 1 z 2 )f2 = 0. (7.13) 

In the variables from ( [7.11| ), equation ( |7,13| ) has form ( |5,22 ) with fj{r) = x(t) — 2. 

In the case \(t) = 8C(x{t) - l)7](t) J T](t)(x(t) - tydt (C = const), particular solutions of (fT.lCQ 
are functions 



w 1 = C{( V (t)) 2 zl - 4z 2 V (t) J r,(t)(x(t) - 3)dt 
w 2 = f(z 1 ,z 2 )exp(l( 1 C) 1 / 2 V (t)z 2 + at)), 

where = -(7C) 1 / 2 / r](t)(x(t) - S)dt + 47C J rj(t)(x(t) - 3)dt)dt and / is an arbitrary 

solution of the following equation 

fi ~ f22 + «X " 2)^- X - (l^r?- 1 + 2( 7 Cy/ 2 )z 2 )f 2 = 0. (7.14) 
After the change of the independent variables 

r = f\ V (t)\exp(4^C) 1 / 2 f V (t)dt)dt, z = \ V {t)\ 1/2 exp(2( 7 C) 1 /2 J v (t)dt)z 2 
in (|7.14|) , we obtain equation (|5.22|) with t){t) = x(t) — 2 again. 



Let us continue to system (7JS). The last equation of (|7.8|) integrates with respect to z 2 to the 
following expression: w 3 = — r]tr]~ 1 z 2 +x- Here x = x(t) is an differentiable function of z\ = 2/3 = t. 
Let us make the transformation from the symmetry group of (|7.2|): 



where ■ $ - f ■ = and 

6 • # + x + mv-Ht-8) - l^r 2 (C • • 0) + M -2 (£ • 0)(0t • 0) = 0. 



Hereafter |^| 2 = $ • This transformation does not modify ansatz (|7.6| ), but it makes 

the function x vanish, i.e., w 3 = —r]t'i]~ 1 z 2 . Therefore, without loss of generality we may assume, 

at once, that w 3 = —r]trj~ 1 z 2 . 

Substituting the expression for w 3 in the other equations of ( [7.8|) , we obtain that 



s = 4H\ 2 (G 2 ^tf^ -ipt-ipt- (ipt ■ ip)rjtv ^l^l 2 + |w7 1 -(vt) 2 v 2 )- 

2^ t -9M- 2 Jw 1 (z 1 ,z 2 )dz 2l 
w{ — r\\i]' 1 z 2 vj\ — |V>| 2 u>22 = 0, 

w\ - r]x'q- 1 z 2 wl - |^| 2 u>| 2 + 7l$l~ 2 (20& ' G)\i>\' 2 z 2 + w^w 2 = 0. 
The change of the independent variables 
t = J(v(tM) 2 dt, z = V (t)z 2 



(7.15) 
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reduces system ( 7,15 ) to the following form: 
w\ ~ w\z = 0> 

, , , (7-16) 

w 2 . - w 2 z + 4 fj 2 (2(ip t ■ 9)fjz + w 1 )w 2 = 0, 

where tJ>(t) = 6{t) = 6(t), t)(t) = r)(t). 

Particular solutions of ( |7.15| ) are the functions 

w l =d + C 2 r,{t)z 2 + C 3 (±( V (t)z 2 ) 2 + J( V (tM) 2 dt), 

w 2 = at, Z2 ) ex P (e(t)z 2 2 + e(t)z 2 + em 

where (£, 2 (t) , S; 1 (t) , £ (t)) is a particular solution of the system of ODEs: 

e t - 2 mr i e - m 2 ^ 2 ? + ^wir 2 = o, 

il - mi' 1 ? - 4^1 W + 27(^ • 0)\M~ A + C 2l rj\i;\- 2 = 0, 

e _ 2 |^|2 e 2 _ |^|2 (C 1 )2 + j{Ci + Cs f(T,(t)$\) 2 dt)$\- 2 = 0, 

and / is an arbitrary solution of the following equation 

fi - \Wf22 + ((Vtr]- 1 + 4H 2 £ 2 )^ + 2\^\ 2 e)f2 = 0. (7.17) 



Equation ( 7.17 ) is reduced by means of a local transformation of the independent variables to the 
heat equation. 

Consider the Lie reductions of system ( |7.2p to systems of ODEs. The second basis operator of 
the each algebra B|, k = 1, 5 induces, for the reduced system obtained from system (|7.2|) by means 
of the first basis operator, either a Lie symmetry operator from Table 2 or a operator giving a 
ansatz of form (|6.4| ). Therefore, the Lie reduction of system ( [7.2[) with the algebras B 2 — B 2 gives 
only solutions that can be constructed for system ( |7.2| ) by means of reducing with the algebras B\ 
and B\ to system (6.1). 

With the algebra B 2 we obtain an ansatz and a reduced system of the following forms: 

v = 0+ A -1 ^ -y)tpi, v 3 = (/> 3 exp(7A((9 1 • y)), 



s = h-\\-\i)\ t .y){9 i -y), 



(7.18) 



(7.19) 



2' 

where $ a = (f) a (uj), h = h(uj), u = t, A = ^V 2 - ip 12 4> 21 = ¥ ■ 9 1 = 4> 2 ■ 2 , 1 = (^ 22 , -4> 21 ) 
P = (-^2^11), and 

4> t + a- 1 ^* • 4>)$\ = °. $ + (i^He 1 ■ 4>) - i 2 \~ 2 {o l ■ e 1 ))^ = o, 

A" 1 ^ • Vt) +Vtr]~ 1 = 0. 
Let us make the transformation from the symmetry group of system (f7^) : 

v(t,y) = v(t,y-i)+it, v 3 (t,y) = v 3 (t,y-£), s(t,y) = s(t,y - f) - ft* • y, 

where 

£ t + \-\e i -o$ + t = Q. (7.20) 

It follows from (7.20) that f tt = \~ l ■ i.e., 0|f£-0'-&t = 0. Therefore, this trasformation 

does not modify ansatz ( 7.18|) , but it makes the functions (j) 1 vanish. And without loss of generality 
we may assume, at once, that (j) 1 = 0. Then 

3 = C , exp(/( 7 A- 1 |0|) 2 dt), C = const. 

The last equation of system ( [7,19| ) is the compatibility condition of system ( |7.2|) and ansatz ( 7.18| ). 
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8 Conclusion 



In this article we reduced the NSEs to systems of PDEs in three and two independent variables 
and systems of ODEs by means of the Lie method. Then, we investigated symmetry properties 
of the reduced systems of PDEs and made Lie reductions of systems which admitted non-trivial 
symmetry operators, i.e., operators that are not induced by operators from A(NS). Some of the 
systems in two independent variables were reduced to linear systems of either two one-dimensional 
heat equations or two translational equations. We also managed to find exact solutions for most of 
the reduced systems of ODEs. 

Now, let us give some remaining problems. Firstly, we failed, for the present, to describe the 



non-Lie ansatzes of form L£ that reduce the NSEs. (These ansatzes include, for example, the well- 
known ansatzes for the Karman swirling flows (see bibliography in jR]]). One can also consider 
non-local ansatzes for the Navier-Stokes field, i.e., ansatzes containing derivatives of new unknown 
functions. 

Second problem is to study non-Lie (i.e., non-local, conditional, and Q-conditional) symmetries 
of the NSEs @. 

And finally, it would be interesting to investigate compatability and to construct exact solutions 
of overdetermined systems that are obtained from the NSEs by means of different additional con- 
ditions. Usually one uses the condition where the nonlinearity has a simple form, for example, the 
potential form (see review |36]]), i.e., rot((«- V)u) = (the NS fields satisfying this condition is 
called the generalized Beltrami flows). We managed to describe the general solution of the NSEs 
with the additional condition where the convective terms vanish |2^, |30fj . But one can give other 
conditions, for example, 

Am = 0, u t + (u- V)n = 0, 

and so on. 

We will consider the problems above elsewhere. 

Appendix 

A Inequivalent one-, two-, and three-dimensional subal- 
gebras of A(NS) 

To find complete sets of inequivalent subalgebras of A(NS), we use the method given, for example, 
in [^7], |28fl . Let us describe it briefly. 

1. We find the commutation relations between the basis elements of A(NS). 

2. For arbitrary basis elements V, W° of A(NS) and each e 6 R we calculate the adjoint action 

W{e) = Ad(eV)W° = Ad(exp(eV))W° 

of the element exp(eV) from the one-parameter group generated by the operator V on W°. This 
calculation can be made in two ways: either by means of summing the Lie series 

oo n 2 

W ( e ) = E -A vn > W °} = W ° + 77 M W °] + "nT [V, [V, W ]} + ..., (A.l) 
z— 1 n\ 1! 21 

n=0 

where {V°, W } = W°, {V n ,W } = [V, {V n ~ l , W }], or directly by means of solving the initial 
value problem 

= \V,W(e)], W(0) = W°. (A.2) 

de 
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3. We take a subalgebra of a general form with a fixed dimension. Taking into account that 
the subalgebra is closed under the Lie bracket, we try to simplify it by means of adjoint actions as 
much as possible. 



A.l The commutation relations and the adjoint representation of the algebra 

A(NS) 

Basis elements (|1.2|) of A(NS) satisfy the following commutation relations: 
[J12, ^23] = — <^3i) [J23, J31] = —J12, [J31, J12] = — J23, 

[d t ,j ab } = [D,j ab } = o, [d t ,D] = 2d t , 

[d t , R(m)] = R(m t ), [D, R(m)] = R(2tm t - m), 

(A.3) 

[dt, Z(x)} = Z(xt), [D, Z(x)} = Z(2t Xt + 2x), 
[R(m),R(n)] = Z(fh tt ■ n - rh ■ n tt ), [J a b, R(m)] = R(rh), 
[J ab ,Zix)} = [Zix),R(m)\ = [Zix),Z( V )\ = 0, 

where m a = m b , m b = —m a , m c = 0, a 7^ b 7^ c 7^ a. 



Note A.l Relations (A.S) imply that the set of operators generates an algebra when, for ex- 
ample, the parameter-functions m a and x belong to C°°((to,ti),M) (Co°((to,ti),W), A((to,ti),R)), 
i.e., the set of infinite- differentiable (infinite- differ entiable finite, real analytic) functions from 
(to,t\) in R, where —00 < to < t\ < +00. But the NSEs admit operators and (1.4 ) 

with parameter-functions of a less degree of smoothness. Moreover, the minimal degree of their 
smoothness depends on the smoothness that is demanded for the solutions of the NSEs Thus, 
ifu a G C 2 ((to,ti) xQ, R) andp £ C 1 ((to 5 ^1) x K) ; where Q is a domain in R 3 , then it is sufficient 
that m a G C 3 ((io, h), R) and x £ C 1 ((to, R)- Therefore, one can consider the "pseudoalgebra" 
generated by operators (jl.^j ). The prefix "pseudo-" means that in this set of operators the commu- 
tation operation is not determined for all pairs of its elements, and the algebra axioms are satisfied 
only by elements, where they are defined. It is better to indicate the functional classes that are 
sets of values for the parameters m a and x in- the notation of the algebra A(NS). But below, for 
simplicity, we fix these classes, taking m a ,x G C°°((to, t\), R) ; and keep the notation of the algebra 
generated by operators (Li) in the form A(NS). However, all calculations will be made in such a 
way that they can be translated for the case of a less degree of smoothness. 
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Most of the adjoint actions are calculated simply as sums of their Lie series. Thus, 



Ad{ed t )D = D + 2ed t , Ad(eD)d t = e~ 2e d t 



Ad{eZ{ x ))d t = d t - eZ( Xt ), Ad(eZ( x ))D = D — eZ(2t X t + 2 X ), 
Ad(eR(m))d t = d t - eR(fh t ) - \e 2 Z(m t ■rh tt -fh- m t tt), 
Ad{eR(m))D = D - eR(2tm t - m)- 

\e 2 Z{2tm t ■ rh u - 2tm ■ m m - 4m • m tt ), 

Ad{eR{m))J ab = J ab - eR(m) + e 2 Z(m a m b tt - m% t m b ), 

Ad(e R(fn))R(n) = R(n) + eZ(m t t ■ n — m ■ n tt ), Ad(sJ a b)R(rn) = R{rh) 

Ad(e Jab) Jed = Jed cos e+ [J ah , J cd ] sine {(a,b) / (c,d) / (b, a)), 
where 

m a = m b , m b = —m a , m c = 0, a / b / c / a, 

fh d = m d cose + fh d sine, m c = m c , a^b^c^a, d£{a;b}. 



(A.4) 



Four adjoint actions are better found by means of integrating a system of form ( A.2| ). As a result 
we obtain that 

M(edt)Z(x(t)) =Z( x (t + e)), Ad(sD)Z( x (t)) =Z(e 2e X (te 2£ )), 

(A.5) 

Ad(edt)R(m(t)) = R(m(t + e)), Ad(eD)R{m{t)) = R(e- e m{te 2e )). 
Cases where adjoint actions coincide with the identical mapping are omitted. 

Note A. 2 // Z(x(t))eA(NS)[C oo ((t ,t 1 ),R)} with -co < t ort\ < +oo, the adjoint representa- 
tion Ad{edt) (Ad(sD)) gives an equivalence relation between the operators Z{x{t)) and Z( X (t + e)) 
(Z( X (t)) and Z(e 2£ x(te 2e ))) that belong to the different algebras 

A(NS)[C°°{{toM)^)] and A(N S)[C°° ({t - e,h - e),R)] 

(A(NS)[C°°((t Q ,ti),R)] and A(NS)[C°°({t e- 2E ,t 1 e' 2e ),R)}) 

respectively. An analogous statement is true for the operator R(m). Equivalence of subalgebras in 



Theorems A.l and A.i is also meant in this sense. 



Note A. 3 Besides the adjoint representations of operators (l.i) we make use of discrete transfor- 



mation (l.t) for classifying the subalgebras of A(NS), 

To prove the theorem of this section, the following obvious lemma is used. 

Lemma A.l Let JV£N. 

A. If X eC N ((to,h),R), then 3 V G C N ((t ,h), R); 2% + 2 V = X . 

B. // X eC iV ((to,ti),M), then 3rj G C N ((t , h),R): 2trj t -r, = X . 

C. If nd G C N ((to,ti),M) andaeR, then 3Z* G C N ((t , h), K): 
2tl\ - I 1 + al 2 = m 1 , 2tl 2 - I 2 - al 1 = m 2 . 
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A. 2 One-dimensional subalgebras 

Theorem A.l A complete set of A(N S)-inequivalent one- dimensional subalgebras of A(NS) is 
exhausted by the following algebras: 

1. A\(k) =< D + 2xJi2 >, where x>0. 

2. A\(h) =< d t + xJ\2 >, where x £ {0; 1}. 

3. Al(n,x) =< Ji2 + R{0, 0,rj(t)) + Z(x(t)) > with smooth functions r\ andx- Algebras A\{rj,x) 
and Al(fj,x) are equivalent if Be, 5 el, BX G C°°((to,ti),R) : 

77(f) = e- e rj(t), x(i) = e 2£ (x(t) + \tt(t)v(t)-\(t)m(t)), (A.6) 
where t = te~ 2e + 5. 

4- A\(rh,x) =< R(m(t)) + Z{x{t)) > with smooth functions m and x'- (^x) ^ (0,0). Algebras 
A{(in,x) and A\(m,x) are equivalent if Be, 5 E R, BC ^ 0, BB £ 0(3), 3fe C°°((to,ii),R 3 )- 

= Ce- £ Bm{t), x(t) = Ce 2£ { X {t) + 4(f) ■ m(t) - m tt (t) ■ f(f)), (A.7) 
where i = te~ 2e + 5. 

Proof Consider an arbitrary one-dimensional subalgebra generated by 
V = aiD + a 2 d t + a 3 J 12 + a 4 J 23 + a 5 J 3 i + i?(m) + 

The coefficients a 4 and 05 are omitted below since they always can be made to vanish by means of 
the adjoint representations Ad(eiJ\ 2 ) and Ad(e 2 J 3 i). 

If ai 7^ we get ai = 1 by means of a change of basis. Next, step-by-step we make a 2 , rh, and 
X vanish by means of the adjoint representations Ad(— ^a 2 a^ 1 dt), Ad(R(i)), and Ad(Z(x)), where 

feC°°(Oo + \a 2 a^i h + ia 2 ar 1 )> R3 )> ??GC°°((to + \a 2 a^ 1 ,fi + i^a^ 1 ), K), 

and I, 7] are solutions of the equations 

2tl t - 1 + a 3 a^(Z 2 , -I 1 , 0) T = fo, 2trj t + 2rj = x + |(4 ■ m - f mtt) 

with m{t) = a^[ l m{t — \a 2 a~{ 1 ) and x(f) = a^xit — \a 2 a^[ l ). Such / and r\ exist in virtue of Lemma 
A.l. As a result we obtain the algebra A\(x), where 2x = a^a^ 1 . In case >t < additionally one 



has to apply transformation ( |1.6| ) with 6 = 1. 

If ai = and a 2 7^ 0, we make a 2 = 1 by means of a change of basis. Next, step-by-step 
we make m and x vanish by means of the adjoint representations Ad(R(l)) and Ad(Z(x)), where 
Te C°°((t ,ti),M 3 ), V e C°°((f ,fi),R), and 

+ G3(^ 2 , -i 1 , 0) T = m, a 2 ??t = x + \{kt -fh-l- m tt ). 

If a 3 = we obtain the algebra ^4 2 (0) at once. If a 3 7^ 0, using the adjoint representation Ad(eD) 
and transformation ( |1.6[) (in case of need), we obtain the algebra ^4 2 (1). 

If a\ = a 2 = and a 3 7^ 0, after a change of basis and applying the adjoint rep- 
resentation Ad(i?(— a% m 2 , a% m 1 , 0)) we get the algebra A^(rj, x)i where r\ = a% 1 m 3 and 
X = a% l X + a 3 2 ( m tt m2 ~ m ^ m tt)- Equivalence relation ( A.6[ ) is generated by the adjoint rep- 
resentations Ad(eD), Ad{5d t ), and Ad(R(0, 0, A)). 

If ai = a 2 = a 3 = 0, at once we get the algebra A\[m, x)- Equivalence relation ( |A.7| ) is generated 
by the adjoint representations Ad(eD), Ad(5dt), Ad(R(l)), and Ad(e a bJ a b)- 
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A. 3 Two-dimensional subalgebras 

Theorem A. 2 A complete set of A(NS)-inequivalent two-dimensional subalgebras of A(NS) is 
exhausted by the following algebras: 

1. A\(x) =< d t , D + xJi 2 >, where x > 0. 

2. A\{k,£) =< D, Ji 2 + i?(0,0,>f|t| 1 / 2 ) + Z{et~ l ) >, where x > 0, e > 0. 

3. Al(x, e) =< d t , Ju + R(0, 0, x) + Z(e) > , where x G {0; 1} ; e > if x=l and 

e G {0;1} i/ x = 0. 

4. Al((T,x,fj,,u,£) =< D + 2xJi2, J R(|t|' T+1/2 (^cosr, z/sinr, //)) + Z(£:|i| CT_1 ) >, w/iere 

r = xln x > 0, /x > 0, v > 0, /i 2 + z^ 2 = 1, ecr = 0, and e > 0. 

5. Al(a, e) =< D, R(0, 0, \t\ a+l / 2 ) + Z(e\t\ a - r ) >, w/iere ecr = and e > 0. 

£. Al(a,fj,,v,e) =< dt + Ju, R{ye at cos t, z^e CT * sint, pe at ) + Z(ee at ) >, where p>0, 
v > 0, p 2 + v 2 = 1, ecr = 0, and £ > 0. 

7. A 2 7 {a, e) =< d t , R(0, 0, e ut ) + ^(ee' 7 *) >, where a G {-1; 0; 1}, ecr = 0, and e > 0. 

5. ^(A,^ 1 ,?,^ 2 ) =< J12 + 12(0,0, A) + Z(^ 1 ), #(0,0,p) + Z(V> 2 ) > with smooth functions (oft) 
X, p, andip 1 : (p, il) 2 ) ^ (0, 0) and X tt p — Xp tt = 0. Algebras A^(X, ip 1 , p, tp 2 ) and A|(A, ip 1 , p, -ip 1 ) 
are equivalent if 3d / 0, 3e, 5, C 2 Gl, 30 G C°°((t ,ii),K): 

X(t) = e £ (X(t) + C 2P (t)), p(t) = C ie - £ p(t), 

ft(t) = e*W\t) + 0tt(t)X(t) - 6(t)X tt (t)+ 

+C 2 (ft(t) + 9 tt (t)p(t)-9(t)p tt (t))), 

ft(t) = C ie 2s (ft(t) + 9 tt (t)p(t) - 9(t) Ptt (t)), 

where t = te~ 2e + 5. 

9. ^(m 1 ,* 1 ,™ 2 ,* 2 ) =< R{m}(t)) + Z(x 1 {t)) 1 R(m 2 (t)) + Z(x 2 (t)) > with smooth functions 
m 1 and x % '■ 

m\ t -m 2 -rh 1 ■ rh 2 tt = 0, rank((m\ x 1 ), {m 2 , x 2 )) = 2. 

Algebras A^{m l , x 1 , m 2 , X 2 ) an d Ag^m 1 , x 1 , rh 2 , x 2 ) are equivalent if 3e,5€M., 3{aij}ij=i >2 ■ 
det{aij}^0, 3BeO(3), 31 £ C°°{(t , h), M 3 ): 

rh l {t) = e~ e aijBfn> '(t), 

(A.9) 

X l {t) = e 2£ aij ( X j (t) + ltt(t) ■ m J '(i) - f(t) ■ n^(t)), 

where t = te~ 2e + 5. 

10. Aj (K, a) =< D + xJi 2 , Zdtj' 7 ) >, w/iere x > 0, <r G R. 
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11. ^(cx) =< d t + Ji 2 , Z(e at ) >, w/iere cr £ K. 

12. A 2 l2 {a) =< d t , Z(e at ) > where a € {-1;0;1}. 



The proof of Theorem A. 2 is analogous to that of Theorem A.l. Let us take an arbitrary 



two-dimensional subalgebra generated by two linearly independent operators of the form 

V i = a\D + aid t + a| J 12 + a|j 23 + a\J 31 + R(rti) + Z( X % 

where a l n = const (n = 1, 5) and [V 1 , V 2 ] £< V 1 , V 2 > . Considering the different possible cases we 
try to simplify V 1 by means of adjoint representation as much as possible. Here we do not present 
the proof of Theorem [A.2| as it is too cumbersome. 



A. 4 Three-dimensional subalgebras 

We also constructed a complete set of J 4(A r 5)-inequivalent three-dimensional subalgebras. It con- 
tains 52 classes of algebras. By means of 22 classes from this set one can obtain ansatzes of 
codimension three for the Navier-Stokes field. Here we only give 8 superclasses that arise from 
unification of some of these classes: 

1. A\ =< D, 8 t , J 12 >. 

2. A\ =< D + xJia, dt, R{0, 0, 1) >, where x > 0. 
Here and below h, a, E\, e 2 , n, v, and a^- are real constants. 

3. A\(o,v,E 1 ,e 2 ) =< D, J12 + v(R{0,0,\t\y 2 ln\t\) + ^(£ 2 |*| _1 In + Z(ei\t\- r ), 
R(0, 0, \t\ a+l / 2 ) + Z{e 2 \t\ a - 1 ) >, where va = 0, e x > 0, v > 0, and ae 2 = 0. 

4. AK^i/^i^a) =< Ji2 + ^(ei) + u(R(0,0,t) + Z{e 2 t)), R{0, 0, e* 7 *) + Z(e 2 e at ) >, 

where va = 0, ae 2 = 0, and, if a = 0, the constants z/, £1, and £ 2 satisfy one of the following 
conditions: 

v = 1, £1 > 0; 1/ = 0, ei = 1, e 2 > 0; ^ = £1 = 0, £ 2 G {0; 1}. 

5. Af^m 1 ,™ 2 ,* 1 ,* 2 ) =<D + 2xJ 12 , R{m l ) + Ztf), R(m 2 ) + Z( X 2 ) >, 
where x> 0, rank(m ,m 2 ) = 2, 

tmj — im l + x(m l2 , —m ll ,0) T = a^m? , 



tXt + X % = aijX 3 , aij = const, 



(an + a22)(a2ifn} ■ m 1 + (a 22 - aujm 1 • m 2 - ai 2 m 2 • m 2 + 

(A 10) 

+2x(m 12 m 21 - m u m 22 )) = 0. 

This superclass contains eight inequivalent classes of subalgebras that can be obtained from it by 
means of a change of basis and the adjoint actions 

Ad(6iD), Ad(5 2 Ji 2 ), Ad(R(n) + Z( V )) 
(Ad(5D), Ad(e ab J ab ), Ad(R(n) + Z(rj))) 
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if a > (x = 0) respectively. Here the functions n and n satisfy the following equations: 
tftt — \n + x(n 2 , — n 1 , 0) T = 6jm*, 

trjt+V = hXi + \t{nttt ■ n - n tt ■ n t ) + n tt ■ n + >i{n x n 2 t - n\ t n 2 ). 

6. Al(K,rn\rh 2 ,x\x 2 ) =<d t + *J 12 , Rim 1 ) + Z{ X V ), R(rh 2 ) + Z( X 2 ) >, 
where ^g{0;1}, rank(m , m 2 ) = 2, 

m t - x(m , -m ,0) = a^m- 7 , t^t = <% x , 



and dij are constants satisfying (A. 10). This superclass contains eight inequivalent classes of sub- 



algebras that can be obtained from it by means of a change of basis and the adjoint actions 
Ad(6id t ), Ad(5 2 J 12 ), Ad(i?(n) + Z(rj)) 

{Adidtdt), Ad(5 2 D), Ad(e ab J ab ), Ad(R(n) + Z{ V ))) 

if x = 1 (x = 0) respectively. Here the functions n and n satisfy the following equations: 
n\ + x(n 2 , — n 1 , 0) T = bifh 1 , 

Vt = biXi + \{n m -n-nu- nt) + *{n x n 2 tt - n\ t n 2 ). 

7. A^(r ] 1 ,r] 2 ,r] 3 ,x) =< J\2 + i?(0, 0, r/ 3 ), R{r] 1 , rj 2 , 0) , i?(-r? 2 , n 1 , 0) >, where 

»7 a eC°°((to,ti),R), r,W-r]Wu = ^ VW^O, 1^0. 
Algebras A 3 (t/ , r/ 2 , r/ 3 ) and y4y (77 1 , ?) 2 , 77 s ) are equivalent if 3S a £ K, 3<54 7^ 0: 

j)l (t) = Sttfit) cos <5 3 - if it) sm5 3 ), 

= ( 5 4 ( r? l(t) s in«53 + 7 ? 2 (t)cOS(53), (A.ll) 

fj 3 (t) = e~ 5l r] 3 (t), 

where t = te~ 2Sl + 5 2 . 

8. A^fh 1 , fn? , m 3 ) =< R{m l ), R(rh 2 ), R(in 3 ) >, where 

fh a e C°°((t ,ii),K 3 ), rank(m 1 ,m 2 ,m 3 ) = 3, m" f • m fe - m a ■ m\ t = 0. 

-* 1 -* 2 -» 3 

Algebras A|(?n , m 2 , m 3 ) and Ag(rh , m ,m ) are equivalent if 35i 6 M 3 , 3B G 0(3), 
3{d a fe} : det{<i a fc} / such that 

m«(f) = d ab Bm b (t), (A.12) 

where i = ie~ 2<51 + 52 . 
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B On construction of ansatzes for the Navier-Stokes field by 
means of the Lie method 

The general method for constructing a complete set of inequivalent Lie ansatzes of a system of 
PDEs are well known and described, for examle, in [^7], |28|]. However, in some cases when the 
symmetry operators of the system have a special form, this method can be modified |]]. Thus, in 
the case of the NSEs, coefficients of an arbitrary operator 

Q = fd t + i a d a + r] a d u a + rfd p 

from A{NS) satisfy the following conditions: 

f0_£0( tj £\ ^o = ^a^ j ^ j n a = rj ab (t,x)u b + ri a0 (t,x), 

(B.l) 

rf = rj 01 (t,x)p + r) 00 (t,x). 

(The coefficients £ a , £°, r] a , and rf also satisfy stronger conditions than (|B.l|) , For example if 
Q G A(NS), then £° = rj ab = const , and so on. But conditions (|B.1| ) are sufficient to 

simplify the general method.) Therefore, ansatzes for the Navier-Stokes field can be constructing 
in the following way: 

1. We fix a M-dimensional subalgebra of A(NS) with the basis elements 

Q m = C°d t + C a d a + (r) mab u b + r) ma0 )d u « + (r] m01 p + r) m00 )d p , (B.2) 

where M e {1; 2; 3}, m = 1,M, and 

rank{(£ m0 ,£ ml ,£ m2 ,£ m3 ), m = TjW} = M. (B.3) 

To construct a complete set of inequivalent Lie ansatzes of codimension M for the Navier-Stokes 
field, we have to use the set of M-dimensional subalgebras from Sec. |A[ Condition ( |B.3| ) is neeeded 
for the existance of ansatzes connected with this subalgebra. 

2. We find the invariant independent variables uj n = u> n (t, x), n = 1, N, where = 4 — M, as 
a set of functionally independent solutions of the following system: 

L m uj = Q m oj = Z m0 d t u> + i ma d a uj = 0, m = T7M, (B.4) 

where L m := t m0 d t + £, ma d a . 

3. We present the Navier-Stokes field in the form: 

u a = f ab (t, x)v b (u>) + g a (t, x), p = f°(t, x)q[p) + g°(t, x), (B.5) 

where v a and q are new unknown functions of <D = {uj n , n = 1, N}. Acting on representation ( |B.5| ) 
with the operators Q m , we obtain the following equations on functions f ab , g a , f°, and g°: 

t m fab m mac fcb Tm„a ^raab J> i „mo0 „ i o 

Ll J —V J i ^ 9 — V 9 -r V ) c— J-j-J) 

(B.6) 

Tin f0 _.m,01 f0 rmJ) „m01 „0 _i_ „m00 

J — V I i ^ 9 — V 9 +V 

If the set of functions f ab , f°, g a , and g° is a particular solution of (|B,6j ) and satisfies the conditions 
rank{(/ lfc , / 2,fc , f 3b ), b = 1,3} = 3 and /° 7^ 0, formulas ( B.5|) give an ansatz for the Navier-Stokes 
field. 
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The ansatz connected with the fixed subalgebra is not determined in an unique manner. Thus, 

if 



{duo 1 
~ — \ + o, 
OUJ n J I n =l,N 



f ab (t, x) = f ac (t, x)F cb {u>), g a (t, x) = g a (t, x) + f ac (t, x)G c {u>), (R7) 
f°(t,x) = f°(t,x)F°(u), ~g°(t,x) = g°(t,x) + f°(t,x)G°(u), 
the formulas 

u a = f a \t, x)v b (ti) + ~g a (t, x), p = f°(t, x)q(0) + ~g°(t, x) (B.8) 



give an ansatz which is equivalent to ansatz (|B.5| ). The reduced system of PDEs on the functions 
v a and q is obtained from the system on v a and q by means of a local transformation. Our problem 
is to find or "to guess", at once, such an ansatz that the corresponding reduced system has a simple 
and convenient form for our investigation. Otherwise, we can obtain a very complicated reduced 
system which will be not convenient for investigation and we can not simplify it. 
Consider a simple example. 

Let M = 1 and let us give the algebra < dt + xJn >, where x £ {0; 1} . For this algebra, 
the invariant independent variables y a = y a (t,x) are functionally independent solutions of the 
equation Ly = (see (B.4)), where 

L := d t + x{x\d X2 - x 2 d Xl ). (B.9) 

There exists an infinite set of choices for the variables y a . For example, we can give the following 
expressions for y a : 

yi = arctan — - xt, y 2 = {x\ + xl) 1/2 , y 3 = x 3 . 

However choosing y a in such a way, for x ^ we obtain a reduced system which strongly differs 
from the "natural" reduced system for x = (the NSEs for steady flows of a viscous fluid in 
Cartesian coordinates). It is better to choose the following variables y a : 

yi = xi cos xt + X2 sin xt, yi = —x\ sin xt + xi cos xt, ys = X3. 

The vector-functions f b = (f lb , f 2b , / 3b ), b = 1, 3, should be linearly independent solutions of the 
system 

Lf l = -xf\ Lf 2 = xf\ Lf = 

and the function /° should satisfy the equation Lf° = and the condition f° 7^ 0. Here the 
operator L is defined by flB,9| ). We give the following values of these functions: 



f 1 = (cosxt,smxt,0), f 2 = (-smxt,cosxt,0), f 3 = (0,0,1), f = 1. 
The functions g a and g° are solutions of the equations 

Lg 1 = -xg 2 , Lg 2 = xg\ Lg 3 = 0, Lg° = 0. 
We can make, for example, g a and g° vanish. Then the corresponding ansatz has the form: 

u 1 = v 1 cos xt — v 2 sin xt, u 2 = v 1 sinxt + v 2 cos xt, u 3 = v 3 , p = q, (B.10) 
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where v a = v a {y\, 1/2, 2/3) and q = q(yi,y2, 2/3) are the new unknown functions. Substituting ansatz 
(B.10) into the NSEs, we obtain the following reduced system: 

qi + xy 2 v\ - xy\v\ - xv 2 = 0, 



v a v\ - vl a 



v a vl 



v a vf, 



r,2 



r,3 



+ <h + xy2V\ - xyiv'i + xv 1 = 



+ <?3 + xy 2 v\ - xyivi = 



r,3 



(B.ll) 



0. 



Here subscripts 1,2, and 3 of functions in ( B.ll ) denote differentiation with respect to y\, y2, and 
2/3 accordingly. System ( |B.11| ), having variable coefficients, can be simplified by means of the local 
transformation 



xy 2 , v = v + xyi, 



(B.12) 



Ansatz ( |B.10| ) and system (B.ll) are transformed under ( B.12| ) into ansatz ( p.2j ) and system ([2.7]). 
where 
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-xx 2 , g = xx!, g 3 = 0, g 



\x\ + xl) 



(B.13) 



71 = — 2x, and 72 = 0. Therefore, we can give the values of g a and g° from (B.13) and obtain 
ansatz (|2.2|) and system (|2.7| ) at once. 

The above is a good example how a reduced system can be simplified by means of modifying 
(complicating) an ansatz corresponding to it. Thus, system (2.7) is simpler than system ( B.ll ) 
and ansatz ( |2.2| ) is more complicated than ansatz ( |B.10| ). 

Finally, let us make several short notes about constructing other ansatzes for the Navier-Stokes 
field. 

Ansatz corresponding to the algebra A\(fh,x) (see Subsec. A.2 ) can be constructed only for 
such t that fh{t) 7^ 0. For these values of t, the parameter-function x can be made to vanish by 
means of equivalence transformations ( |A.7| ). 

Ansatz corresponding to the algebra ^4§(A, ip 1 , p, 1I) 2 ) (see Subsec. A.3 ) can be constructed only 
for such t that p(t) ^ 0. For these values of t, the parameter-function tp 2 can be made to vanish by 
means of equivalence transformations ( |A.8| ). Moreover, it can be considered that Xtp — Xpt £ {0; 1}. 
The algebra obtained finally is denoted by ^4g(A, x, P, 0). 

Ansatz corresponding to the algebra A^m 1 , x 1 , m 2 , X 2 ) ( see Subsec. A3) can be constructed 
only for such t that rank(m 1 ,m 2 ) = 2. For these values of t, the parameter-functions x l c& n be 
made to vanish by means of equivalence transformations ( |A.9f) . 

The algebras A 2 (x,a), A^i^a), and A 2 2( a ) can n °t be used to construct ansatzes by means 
of the Lie algorithm. 

In view of equivalence transformation ( A.ll| ), the functions rf in the algebra A?(r] , rj 2 , rj 3 ) (see 
Subsec. A. 4) can be considered to satisfy the following condition: 
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